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Structure of the Thesis 

The main object studied in this thesis, in the differential Galois theory framework, 
is the one-dimensional stationary non-relativistic Schrodinger equation 

dl^ = (V - A)*, V £ K, 

where K is a differential field containing x, closed algebraically and of characteris- 
tic zero. In particular, K is considered as the smallest differential field containing 
the potential V. 

This thesis is divided in two parts. 
Chapter 1. Theoretical Background 

In this part there are not original results. Summaries of Picard-Vessiot theory and 
supersymmetric quantum mechanics, necessaries to understand the next chapter, 
are presented here. 

Chapter 2. Differential Galois Theory Approach to Supersymmetric 
Quantum Mechanics 

This part contain the original results of this thesis, which were developed using 
the previous chapter. Up to specific cases; theorems, propositions, corollaries and 
lemmas given in this chapter are considered as original results of this thesis. 

Two different Galoisian approaches are studied in this chapter, which de- 
pends on the differential field: the first one is C(x) and the second one is K = 
C(z(x),d z (x)), where z = z(x) is a Hamiltonian change of variable. This concept 
allow us to introduce an useful derivation d z , important tool to transforms dif- 
ferential equations with non rational coefficients into differential equations with 
rational coefficients, to apply the results given in the case of C(x). 
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This chapter is divided in three parts or sections. 



Section 2.1 

Here is introduced the set A as the set of values of A in which the Schrodinger 
equation is integrable in the sense of Picard-Vessiot theory. With this new set A, 
we define in an easy way the concepts of algebraically solvable, algebraically quasi- 
solvable and algebraically non-solvable potentials. 

Also we define in this section the concepts of iso-galoisian, virtually iso- 
galoisian and strong iso-galoisian transformations. So, we give in proposition 2.1.5 
some conditions in which the transformation of the second order linear differential 
equation into a reduced form (the term in d x y is absent) is strong iso-galoisian 
or virtually strong iso-galoisian. In particular case, corollary 2.1.7 shows that the 
reduction of Sturm-Liouville problem is a virtually strong iso-galoisian transfor- 
mation. 

Section 2.2 

This section is devoted to the analysis of supersymmetric quantum mechanics 
from a differential Galois theory point of view when V is a rational function. We 
start considering the polynomial case proving that the only one possibility to get 
algebraically solvable potential is when V is a polynomial of degree 2. In case of 
algebraically quasi-solvable potentials are presented as examples the well known 
quartic and sextic anharmonic oscillators. Some results such as lemma 2.2.1 and 
theorem 2.2.2 were previously published in [5] and used in [?]. 

Kovacic's algorithm is applied to solve the some Schrodinger equations with 
rational potential (shape invariant potentials) obtaining also their differential dif- 
ferential Galois groups and cigenrings. In proposition 2.2.6 and corollary 2.2.7 is 
shown that the Schrodinger equation with rational potential does not fall in case 
3 of Kovacic's algorithm. 

Darboux transformation and Crum iteration are written in a Galoisian sense 
(theorem 2.2.8, proposition 2.2.11 and proposition 2.2.12), so that by proposition 
2.2.9 and proposition 2.2.10, such transformations can be seen as isogaloisian and 
for instance they preserves the Eigenrings. We prove in proposition 2.2.13 that the 
supersymmetric partner potentials and the supcrpotcntial are rational functions 
and in particular, corollary 2.2.14 says that the Darboux transformation is strong 
iso-galoisian when the supcrpotcntial is a rational function. 

Finally, we define, in the context of differential Galois theory, the concept 
of shape invariant potentials giving an algorithm to check whether a potential 
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satisfy the shape invariance condition (remark 2.2.16) and providing an important 
result (theorem 2.2.17) in where the Schrodingcr equations obtained with the 
procedure of shape invariant potentials preserves the differential Galois groups 
and the Eigcnrings. 

Section 2.3 

This section is devoted to the supersymmetric quantum mechanics with non- 
rational potentials. We start considering the proper pullback between two second 
order linear differential equations in which one of them has rational coefficients 
(algebrization process when the first differential equation with non-rational coeffi- 
cients is transformed into a differential equation with rational coefficients). Propo- 
sition 2.3.2 shows that the change of independent variable preserves the connected 
identity component of the differential Galois group and give one formula of the 
transformed equation (see also [5]). 

In proposition 2.3.3, are presented the relationships between the differen- 
tial fields, Picard-Vessiot extensions and differential Galois groups corresponding 
to this change of independent variable. Using these propositions, in remark 2.3.4 
we introduce the concept of Hard algebrization giving an algorithm to obtain it 
whether is possible. Proposition 2.3.5, which also appears in [5], shows the al- 
gebrization of some families of second order differential equations in where the 
differential field is given by C(x, e$ f), being / £ C(x). 

The definition of Hamiltonian change of variable is introduced, leading to 
the concept of Hamiltonian algebrization, which allows to arrive to one systematic 
algebrization of second order linear differential equations presented in proposition 
2.3.9, remark 2.3.10 and proposition 2.3.11 (these results also can be found in 
[5] and were applied in [H[5]). We show through examples that the algebrization 
process can be used as transformation between differential equations without ex- 
pecting to obtain differential equations with rational coefficients. 

We introduce a new derivation d z = ^fud Zl in where z = z(x) is a Hamil- 
tonian change of variable and d x z = y/a. The most important theoretical result 
here is that the differential field K = C(z(x),d x z(x)) is isomorphic to the differ- 
ential field K = C(z, yfa) D C(z), which allows to preserve the differential Galois 
groups and Eigcnrings in the algebrization process (theorem 2.3.13 and propo- 
sition 2.3.16). So, we transform a lot of differential equations with non-rational 
coefficients into differential equations with rational coefficients (Riccati, systems, 
etc..) and in the case of second order linear differential equations, Kovacic's algo- 
rithm can be used successfully. 

We recover the results of section 2.2 rewriting everything with ^ (hat over 
the symbols used in section 2.2), in this way the algebrized supersymmetric quan- 
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turn mechanics is presented with its elements: algebrized Schrddinger operator 
H, algebrized superpotential W, algebrized super symmetric partner potentials V±, 
algebrized shape invariant potentials, algebrized Darboux transformation DT, al- 
gebrized Crum iteration CI„, algebrized ladder operators A, A\ algebrized wave 
functions 4", etc.. An important fact is given in theorem 2.3.20, i.e., = i^DT, 
where ip is the Hamiltonian algebrization. 

Using the Hamiltonian algebrization and Kovacic's algorithm we solve some 
Schrodinger equations with non-rational potentials (shape invariant potentials), 
obtaining the differential Galois groups and the Eigenrings of these differential 
equations. 

Finally, we give a mechanism to search new exactly solvable potentials through 
the inverse process of the Hamiltonian algebrization, using known parameterized 
differential equations. 

Historical Outline 

This historical outline begins with two mathematicians: Gaston Darboux and 
Emilc Picard. Darboux published in 1882 the paper [35] in where he presents 
a proposition in a general way, which in particular case the history proved to be 
a notable theorem today known as Darboux transformation. Darboux had shown 
that whenever one knows to integrate the equation 

d 2 x y = (f(x) + m)y 

for all the values of the constant m, one can obtain an infinite set of equations, 
displaying the variable parameter in the same way. which arc intcgrablc for any 
value of the parameter. This proposition also can be found in his book [29l p. 210]. 

One year after, in 1883, Picard published the paper [73] in which he gave 
the starting point to a Galois theory for linear differential equations. Although 
the analogies between the linear differential equations and the algebraic equations 
for a long time were announced and continued in different directions, Picard de- 
veloped an analogue theory to the Galois theory for algebraic equations, arriving 
to a proposition which seems to correspond to the fundamental Galois theorem, 
in where he introduces the concept of group of linear transformations correspond- 
ing to the linear differential equation, which today is known as Differential Galois 
Group (the group of differential automorphism leaving fixed the elements of the 
field base). Another contribution of Picard to this Galois theory was the paper 
[72] in 1887. 

Five years after, in 1892, Ernest Vessiot, doctoral former student of Picard, 
published his thesis [103j giving consolidation to the new Galois theory for linear 
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differential equations, the so-called Picard-Vessiot theory. Two years later, in 1894, 
Picard published the paper [73] , summarizing the results presented in [721 [751 1103] 
which also can be found in his book [751 §7 ]. 

Curiously, Picard-Vessiot theory and Darboux transformation were forgotten 
during decades. The Picard-Vessiot theory was recovered by Joseph Fels Ritt (in 
1950, see [85]), Irving Kaplansky (in 1957, see [53] ), and fundamentally by Ellis 
Kolchin (in 1948, see [55] and references therein). Kolchin wrote the Differential 
Galois Theory in a modern language (algebraic group theory). 

Darboux transformation was presented as an exercise in 1926 by Ince (see 
exercises 5, 6 and 7 [49l p. 132]), which follows closely the formulation of Darboux 
given in [251 [29] . 

In 1930, P. Dirac publishes The Principles of Quantum Mechanics, in where 
he gave a mathematically rigorous formulation of quantum mechanics. 

In 1938, J. Delsarte wrote the paper [30], in which he introduced the notion 
of transformation (transmutation) operator, today know as intertwining operator 
which is closely related with Darboux transformation and ladder operators. 

In 1941, E. Schrodingcr published the paper [89] in which he factorized in 
several ways the hypergeometric equation. This was a byproduct of his factoriza- 
tion method originating an approach that can be traced back to Dirac's raising 
and lowering operators for the harmonic oscillator. 

Ten years later, in 1951, another factorization method was presented. L. In- 
feld and T. E. Hull published the paper [50] in where they gave the classification 
of their factorizations of linear second order differential equations for eigenvalue 
problems of wave mechanics. 

In 1955, M.M. Crum inspired in the Liouville's work about Sturm-Liouvillc 
systems (see [61] [62]), published the paper [27] giving one kind of iterative gener- 
alization of Darboux transformation. Crum surprisingly did not mention Darboux. 

In 1971, G.A. Natanzon published the paper [69], in which he studied a gen- 
eral form of the transformation that converts the hypergeometric equation to the 
Schrodingcr equation writing down the most general solvable potential, potential 
for which the Schrodingcr equation can be reduced to hypergeometric or confluent 
hypergeometric form, a concept introduced by himself. 

Almost one hundred years later than Darboux's proposition, in 1981, Edward 
Witten in his renowned paper |107j gave birth to the Super symmetric Quantum 
Mechanics, discussing general conditions for dynamical supersymmetry breaking. 
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Since the work of Witten, thousands of papers, about supcrsymmctric quan- 
tum mechanics, has been written. We mention here some relevant papers. 

In 1983, L. E. Gcndcnshtcin published the paper |37j in where the Shape in- 
variance condition, i.e. preserving the shape under Darboux transformation, was 
presented and used to find the complete spectra for a broad class of problems in- 
cluding all known exactly solvable problems of quantum mechanics (bound state 
and reflectionless potentials). Today this kind of exactly solvable potentials satis- 
fying the shape invariance condition are called Shape invariant potentials. 

In 1986, A. Turbiner in |102j introduces the concept of quasi- exactly solvable 
potentials, giving an example that is well known as Turbiner 's potential. 

In 1991, V.B. Matveev and M. Salle published the book in where they 
focused on Darboux transformations and their relation with solitons. Matveev 
and Salle interpreted the Darboux transformation as Darboux covariance of a 
Sturm-Liouvillc problem and also proved that Witten's supersymetric quantum 
mechanics is equivalent to a single Darboux transformation. 

In 1996, C. Bender and G. Dunne studied the sextic anharmonic oscilla- 
tor in [11) . which is a quasi-exactly solvable model derived from the Turbiner's 
potentials. They found that a portion of the spectrum correspond to the roots 
of polynomials in the energy. These polynomials are orthogonal and are called 
Bender-Dunne polynomials. 

Relationships between the spectral theory and differential Galois theory have 
been studied by V. Spiridonov [97], F. Beukers [16] and Braverman et. al. [19]. As 
far as we know, Spiridonov was the first author that considered the usefulness of 
the Picard-Vessiot theory in the context of the quantum mechanics. This thesis 
agrees with his point of view. 



Chapter 1 

Theoretical Background 



In this chapter we set the main theoretical background needed to understand the 
results of this thesis. We start setting conventions and notations that will be used 
along this work. 

• The sets Z + , Z_ , Zl and Zl are defined as 

Z+ = {neZ: n>0}, Z_ = {n E Z : n < 0}, Z* + = Z+, Z*=Z". 

• The cardinality of the set A will be denoted by Card(A). 

• The determinant of the matrix A will be denoted by det A. 

• The set of matrices n x n with entries in C and determinant non-null, the 
general linear group over C, will be denoted by GL(n, C). 

• The derivation d/d£ will be denoted by d^. For example, the derivations 
' = d/dx and ' = d/dt are denoted by d x and dt respectively. 

1.1 Picard-Vessiot theory 

Picard-Vessiot theory is the Galois theory of linear differential equations. In the 
classical Galois theory, the main object is a group of permutations of the roots, 
while in the Picard-Vessiot theory it is a linear algebraic group. For polynomial 
equations we want a solution in terms of radicals, which from classical Galois 
theory it is well if the Galois group is a solvable group. 

An analogous situation holds for linear homogeneous differential equations 
(see p.5, 26, 67, 81 ). The following definition is true in general dimension, but for 
simplicity we are restricting to matrices 2x2. 



8 



Chapter 1. Theoretical Background 



1.1.1 Definitions and Known Results 

Definition 1.1.1. An algebraic group of matrices 2 x 2 is a subgroup G C GL(2, C), 
defined by algebraic equations in its matrix elements and in the inverse of its 
determinant. That is, for A € GL(2,C) given by 

, / in 112 \ , , „ 

A=[ , det A = X\\X22 — X2\X22 

\ X 2 1 X 2 2 J 

there exists a set of polynomials 



{Pi(x U , Xi2,X21,X 2 2, 1/ detA)} ieI , 

such that 



AeG ^ Vi £ I,Pi(xn,xi2,X2i,x 2 2A/ det A) = 0. 



In this case we say that G is an algebraic manifold endowed with a group 
structure. 

Examples (Known algebraic groups). The following algebraic groups should be 
kept in mind throughout this work. 



• Special linear group group 
b 

c d 



SL(2,C) = <j( a b ), ad-bc = l, a,b,c,deC 



I c" 1 i 

( g ^ > d e C 
Infinite dihedral group (also called meta-abelian group): 

• n— quasi-roots: — | c^ 1 ^ ' c n — 1, d £ 

1 c" 1 1 



• Identity group: e 



1 
1 
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• The tetrahedral group A 4 2 of order 24 is generated by matrices 

Mi = and M 2 = 1(2£ - 1) ^ ^ , 

where £ denotes a primitive sixth root of unity, that is, £ 2 — £ + 1 = 0. 

• The octahedral group S^ 2 of order 48 is generated by matrices 

Mi = (o and M 2 = ^ a + i)(J _\), 

where £ denotes a primitive eighth root of unity, that is, £ 4 + 1 = 0. 

• The icosahedral group Af h2 of order 120 is generated by matrices 

Mo c') ™ d -*)' 

where £ denotes a primitive tenth root of unity, that is, £ 4 — £ 3 +£ 2 — £ + 1 = 0, 

= £3 _ £2 +4 £ _ 2 and ^ = £3 + 3^2 _ 2 ^ + L 

Recall that a group G is called solvable if and only if there exists a chain of 
normal subgroups 

e = Go < Gi < . . . < G n = G 

such that the quotient Gi/Gj is abclian for all n > i > j > 0. Also recall that 
an algebraic group G has a unique connected normal algebraic subgroup G° of 
finite index. This means that the identity connected component G° is the largest 
connected algebraic subgroup of G containing the identity. For instance, if G = G°, 
we say that G is a connected group. 

Furthermore, if G° satisfy some property, then we say that G virtually satisfy 
such property. In this way, virtually solvability of G means solvability of G° and 
virtually abelianity of G means abelianity of G° (see |106j ). 

Theorem 1.1.2 (Lie-Kolchin). Let G C GL(2, C) be a virtually solvable group. 
Then G° is triangularizable, that is, conjugate to a subgroup of upper triangular 
matrices. 

Definition 1.1.3. Let G C GL(2, C) be a group acting on a vector space V. We say 
that (the action of) G is either: 

1. Reducible, if there exists a non-trivial subspace W C V such that G(W) C W. 
We say that G is irreducible if G is not reducible. 

2. Imprimitive, if G is irreducible and there exists subspaces Vi such that V — 
V\ ®- ■ •®V m , where G permutes transitively the Vi, i.e Vi = 1, . . . , to, Wg G G, 
3j G {1, . . . , to} such that g{Vi) = Vj. We say that Vi, . . . , V m form a system 
of imprimitivity for G. 
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3. Primitive, if G is irreducible and not imprimitive. 

Examples. Any subgroup of the Borel group is reducible, the infinite dihedral 
group is imprimitive and the groups ^4 SL2 , Sf h2 , Af h2 , SL(2,C) are primitives 
(see |5T1ITP5]L 

Definition 1.1.4 (Differential Fields). Let K (depending on x) be a commutative 
field of characteristic zero, d x a derivation, that is, a map d x : K — > K satisfying 
d x (a + b) = d x a + d x b and d x (ab) = d x a -b + a- d x b for all a, b £ K. By C we denote 
the field of constants of K 

C = {c£ K\d x c = 0} 

which is also of characteristic zero and will be assumed algebraically closed. In 
this terms, we say that K is a differential field with the derivation d x . 

Along this work, up to some specifications, we consider as differential field 
the smallest differential containing the coefficients. Furthermore, up to special 
considerations, we analyze second order linear homogeneous differential equations, 
that is, equations in the form 

C := d 2 x y + ad x y + by = 0, a,b£K. (1.1) 

Definition 1.1.5 (Picard-Vcssiot Extension). Suppose that 2/1,2/2 is a basis of solu- 
tions of C given in equation (jl.lj) . i.e., yi,?/2 arc linearly independent over K and 
every solution is a linear combination over C of these two. Let L = K(y\,y2) = 
K(yi,y2,d x yi, d x y2) the differential extension of K such that C is the field of con- 
stants for K and L. In this terms, we say that L, the smallest differential field 
containing K and {2/1,2/2}, is the Picard-Vessiot extension of K for C. 

Definition 1.1.6 (Differential Galois Groups). Assume K, L and C as in previ- 
ous definition. The group of all differential automorphisms (automorphisms that 
commutes with derivation) of L over K is called the differential Galois group of 
L over K and is denoted by DG&\(L/K). This means that for a £ T)G&\(L/K), 
a(d x a) = d x (o~(a)) for all a £ L and Va £ K, <r(a) = a. 

Assume that {j/i, 1/2} is a fundamental system of solutions (basis of solutions) 
of C. If cr £ DGal(L/-ftT) then {ayi, cry2] is another fundamental system of C. Hence 
there exists a matrix 

K= (° J)eGL(2,C), 

such that 




in a natural way, we can extend to systems: 

\d x yi d x y 2 ) \o{d x yi) a(d x y 2 ) J \d x yi d x y 2 ) 



1.1. Picard-Vcssiot theory 
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This defines a faithful representation DGal(L/K) — > GL(2,C) and it is pos- 
sible to consider DGal(L/_ftT) as a subgroup of GL(2,C). It depends on the choice 
of the fundamental system {2/1,2/2}, but only up to conjugacy. 

One of the fundamental results of the Picard-Vessiot theory is the following 
theorem (see [551 [55] ). 

Theorem 1.1.7. The differential Galois group DGal(L / K) is an algebraic subgroup 
o/GL(2,C). 

Examples. Consider the following differential equations: 

• C := d x y — 0, the basis of solutions is given by yi — 1, 2/2 = x. If we set 
as differential field K = C(x), we can see that <r(l) = 1, o~(x) = x, then the 
Picard-Vessiot extension L = K and for instance DGdl(L/K) = e: 

'(£)-<» < !)-&)■ 

Now, if we set K — C, then L = K(x), d x x G C, d x (a(x)) = cr(d x x) = cr(l) = 
1 = d x x, so a(x) = x + d, d € C and for instance DGal(L/ii") = <G a : 




• £ := d x y = ny, k G C*, the basis of solutions is given by j/i = e^ x , 
2/2 = e~^ x , with k ^ 0. If we set as differential field K = C(x), we can see 
that L = K(e^ x ) = K(e^ x ), 

f d x yA _ d x {a{yi)) _ d x yi f d x y 2 \ _ d x (a(y 2 )) _ d x y 2 

V yi J o'ivi) yi \ j/2 ) o-(y 2 ) 2/2 

o-{yiy<z) = o-(y 1 )o-(y 2 ) = 2/12/2 = 1, cr(yi) = qji, <r{y 2 ) = dy 2 , c, d G C, but 
cd = 1 and for instance DGal(L/A') = G,„: 




Now, if we set K = C, we obtain the same result. 
• C := d x y + I ^-9 X 2/ = 0, the basis of solutions is given by 2/1 = z, where 
z n = x, 2/2 = 1. If we set K = C(x), then L = k(x»\, y{ 1 = x G C(x), 

o~ n (yi) — a (yi) = x, a (yi) = c 2/i> so th & t c" = 1 and for instance DGal(L/AT) 
is given by: 
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• C := d 2 y + 4^2 V = 0, n e Z, the basis of solutions is given by j/i = a; "at 1 1 
j/2 = x J ^' . If we set Jv = C(x) and n even, then L = K (x^\ , 

tr(y 1 ) = cy 1 , o^iVi) = c 2n yf n = cr(y* n ) = yj n , c 2n = 1, 

a(y 2 ) = dy 2) a 2n (y 2 ) = d 2n y 2 2 n = *(y 2n ) = y 2 2 n , d 2n = 1, 

c(yiy2) = 2/12/2 = &{yi)o~{y 2 ) = cdy\y 2 so that cd = 1 and for instance 
DGal(L/iT) = G[ 2n l: 




= 1, n > 1. 



Now, if we consider n odd, then L = K (x«) , and DGal(L/if) = GPH 
• Cauchy-Euler equation 

„o m(m + 1) „ 
C:=a£y= V 2 ^ meC, 

the basis of solutions is y\ = x m+1 , y 2 = x~ m . Setting K = C(x), we have 
the following possible cases: 

- for meZ, L = K and DGal(L//\) = e, 

- for m e Q \ Z, L = K{x m ) and DGal(L/K) = G [d] , where m = n/d, 

- for meC\Q,i = X(:r m ) and DGal(i/A") = G m . 

Definition 1.1.8 (Intcgrability). Consider the linear differential equation C such 
as in equation (jl.ip . We say that C is integrable if the Picard-Vessiot extension 
L D K is obtained as a tower of differential fields K = Lq C £1 C ■ ■ • C L TO = i 
such that L; = Li-iirf) for i = 1, . . . , m, where either 

1. 77 is algebraic over £i_i, that is 77 satisfies a polynomial equation with coef- 
ficients in Li— 1. 

2. 77 is primitive over Lj_i, that is € 

3. 77 is exponential over that is d x rj/rj G £,_i. 

We recall that the differential field of coefficients has been fixed before, i.e., 
the smallest differential field containing the coefficients. 

We remark that the usual terminology in differential algebra for integrable 
equations is that the corresponding Picard-Vessiot extensions are called Liouvil- 
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Theorem 1.1.9 (Kolchin). The equation C given in is integrable if and only 

ifDGal(L/K) is virtually solvable. 

Consider the differential equation 

£:=^C = rC, reK. (1.2) 

We recall that equation fll.2|) can be obtained from equation (jl.ip through 
the change of variable 

¥ = r=^ + ^-6 (1-3) 

and equation (| 1 . 2[) is called the reduced form of equation (jl.ip . 

On the other hand, introducing the change of variable v = d x (/( we get the 
associated Riccati equation to equation (|1.2p 

d x v = r — v 2 , v = (1.4) 

where r is obtained by equation (|1.3[) . 

Theorem 1.1.10 (Singer 1981, [H]). T/ie Riccatti equation (|1.4[) fras one algebraic 
solution over the differential field K if and only if the differential equation fj 1 . 2|) 
is integrable. 

For £ given by equation (|1.2j) . it is very well known (see [53l EH [81]) that 
DGal^(£) is an algebraic subgroup of SL(2, C). The well known classification of 
subgroups of SL(2,C) (see [53l p.31], [HI p.7,27]) is the following. 

Theorem 1.1.11. Let G be an algebraic subgroup o/SL(2,C). Then, up to conju- 
gation, one of the following cases occurs. 

1 GC1 and then G is reducible and triangular izable. 

2. G $4 B, G C and then G is imprimitive. 

3. G G {A^ 1 " 2 , 5*4 L2 , ^5 L2 } and then G is primitive (finite) 

4- G = SL(2,C) and then G is primitive (infinite). 

Definition 1.1.12. Consider the differential equation C given by equation ()1.2|) . 
Let {Ci>C2} be a fundamental system of C. Let / = /(Yi,!^) G C[Yi,i2] be a 
homogeneous polynomial, we say that: 

1. The polynomial / is an invariant with respect to C if its evaluation on a C- 
basis {(1, £2} of solutions is invariant under the action of DGal(i//C), that is, 
for every tr G DGal(i/X), 0-/1(2;) = /i(.t), where /i(.t) = /(Ci(z), C, 2 {x)) G A". 
The function h{x) is called the m/we of the invariant polynomial /. 
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2. The polynomial / is a semi-invariant with respect to C if the logarithmic 
derivative %^ of its evaluation h(x) = f(Ci{%), C2( x )) on any C-basis {Ci, C2} 
is an element of K, that is, for every a G T)Gai(L/K), oO = 9, where 9 = 
d x h(x)/h{x) G K. 

Theorem 1.1.13 (Kovacic, |56|). Let {Ci,C2} & e a fundamental system of solutions 
of C given by the differential equation (|1.2p . Then, for some i G {1,2} and for 
every a G DGal(L / ' K) , exclusively one of the following cases holds. 

1. DGal(L / K) is reducible and then f — Yi is semi-invariant with respect to £, 
i.e d x {\nQ G K. 

2. Y>G'a\(L / K ) is imprimitive and then f\ — is semi-invariant with respect 
to C, f% = (Y1Y2) 2 is invariant with respect to C, i.e [K (d x (In £,•)) : K] = 2. 

3. DGal(L / ' K) is finite primitive and then the invariants with respect to C is 
either h = (Yf + 8^1^) 3 , or f 2 = {Y(?Y 2 - YiY 2 5 ) 2 or f 3 = Y?Y 3 - 
UY?Y$ - Y X Y^, i.e [K (d x (In CO) : K] = 4,6, 12 . 

4- DGal(L/_ftT) is infinite primitive, i.e there are no non-trivial semi-invariants. 

Statements and proofs of theorems 11.1.91 II. 1.101 and II. 1.111 can be found in 

1.1.2 Kovacic 's Algorithm 

Considering K = C(x), C = C in theorems 11.1.91 fl. 1.111 and 11.1.131 Kovacic in 
1986 ([56]) introduced an algorithm to solve the differential equation (| 1 . 2[) show- 
ing that (|1.2j) is integrable if and only if the solution of the Riccati equation (|1.4p 
is a rational function (case 1), is a root of polynomial of degree two (case 2) or is 
a root of polynomial of degree 4, 6, or 12 (case 3). For more details see reference 
|56j . Improvements for this algorithm are given in references |331 1341 1101) . Here, 
we follow the original version given by Kovacic in reference |56j with an adapted 
version presented in reference [5] . 

Each case in Kovacic's algorithm is related with each one of the algebraic 
subgroups of SL(2, C) and the associated Riccatti equation 

d x v — r — v 2 = (yfr — v) (\^r + v) , v = 

According to Theorem 11.1.111 there are four cases in Kovacic's algorithm. 
Only for cases 1, 2 and 3 we can solve the differential equation, but for the case 4 
the differential equation is not integrable. It is possible that Kovacic's algorithm 
can provide us only one solution so that we can obtain the second solution 
(C2) through 

C2 = Ci / J- (1-5) 
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Notations. For the differential equation given by 
^C = rC, r = j, s,teC[x), 

we use the following notations. 

1. Denote by V be the set of (finite) poles of r, V = {c G C : t(c) = 0}. 

2. Denote by V = V U {oo}. 

3. By the order of r at c G T', °(r c ), we mean the multiplicity of c as a pole of 
r. 

4. By the order of r at oo, o (r^) , we mean the order of oo as a zero of r. That 
is o (r^ = deg(i) - deg(s). 

The four cases 

Case 1. In this case [y/r] c and [y^loo means the Laurent series of ypr at c and 
the Laurent series of y/r at oo respectively. Furthermore, we define e(p) as follows: 
if p G r, then e (p) e {+,—}. Finally, the complex numbers a+ , a~ , a+ , will 
be defined in the first step. If the differential equation has no poles it only can fall 
in this case. 

Step 1. Search for each c G V and for oo the corresponding situation as 
follows: 



(c ) If o (r c ) = 0, then 
(ci) If o (r c ) = 1, then 
(C2) If o (r c ) = 2, and 



[VF] c = o, «± = o. 

[VF] c = 0, a± = l. 

• + b(x — c)~ 2 + • • • , then 



a, = 



(c 3 ) If o(r c ) = 2v > 4, and 

r = (a (x - c)~" + ... + d(x- c)~ 2 ) 2 + 6(21 - c)~ {v+1) + ■■■ , then 

[v^] c = a(x- c)~ v + ... + d(x- c)~ 2 , af = i (±- + w] . 
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± 1 ± Vl + 46 



(001) If o (foe) > 2, then 

[\/r] oo =0 ) a+=0, a-=l. 

(002) If (j"oo) = 2, and r = • • • + 6a; 2 + • • • , then 

(003) If o (roo) = — 2v < 0, and 

r = (ax v + ... + df + fa-"" 1 + • • • , then 

[V^loo = a2;t ' + ■■• + d > and a to = \ f±^ - « 

Step 2. Find D^U defined by 

D = jn G Z + : „ = a^°°) - X>c (c) ,V( e • 

If D = 0, then we should start with the case 2. Now, if Card(-D) > 0, then for 
each n £ D we search tu £ C(x) such that 

w = e (00) [V^ + £ (e W [Vr] c + «; (c) (* - c)" 1 ) . 
cer' 

Step 3. For each n £ D, search for a monic polynomial P n of degree n with 

d 2 P„ + 2u>d x P n + {d x oj + oj 2 - r)P n = 0. (1.6) 

If success is achieved then £1 = P n eJ w is a solution of the differential equation. 
Else, case 1 cannot hold. 

Case 2. Search for each c £ T' and for 00 the corresponding situation as 
follows: 

Step 1. Search for each c £ V and 00 the sets E c ^% and E x 7^ 0. For each 
ceT' and for 00 we define E c C Z and C Z as follows: 

(ci) If o (r c ) = 1, then E c = {4} 

(C2) If o (r c ) = 2, and r = ■ ■ ■ + b(x — c)" 2 + • • • , then 



|2 + fcVl + 46:/c = 0,±2| 
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(c 3 ) If o (r c ) = v > 2, then S c = {v} 

(001) If o (roo) > 2, then = {0, 2, 4} 

(002) If i r oo) = 2, and r = ■ ■ ■ + bx 2 + ■ ■ ■ , then 



1 2 + kVTTib :k = 0, ±2 1 . 



(oo 3 ) If o (r^) = v < 2, then = {v} 
Step 2. Find D^U defined by 



^ ^eoo - ^e c ^ , Ve p G £"p, pGrj 



£) = G Z+ : n = 2 

If -D = 0, then we should start the case 3. Now, if Card(D) > 0, then for each 
tiEf we search a rational function defined by 



9 Z— ' i™ 



2 ^ x - c 
cer' 



Step 3. For each n £ D, search a monic polynomial P„ of degree n, such that 

dlP n +WdlP„ + {3d x e+W 2 -Ar)d x P n + (d x 2e + 36d x 6 + (9 3 - 4r<9 - 25 T r) P„ = 0. 

(1.7) 

If P n does not exist, then case 2 cannot hold. If such a polynomial is found, set 
4> = + d x P n /P n and let u> be a solution of 



2 + ^ + 1 (d x <j> + 2 - 2r) = 0. 



2 



Then £1 = e-> u is a solution of the differential equation. 

Case 3. Search for each c G V and for 00 the corresponding situation as 
follows: 

Step 1. Search for each c G V and 00 the sets E c ^% and E^ ^ 0. For each 
ceT and for 00 we define E c C Z and C Z as follows: 

(ci) If o (r c ) = 1, then £ c = {12} 

(c 2 ) If o (r c ) = 2, and r = • • • + 6(2; — c)~ 2 + • • • , then 

£ c = 1 6 + fcVl + 46 : fc = 0, ±1, ±2, ±3, ±4, ±5, ±6} . 
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(oo) If o (roo) = u > 2, and r = • • • + bx 2 + • • • , then 

6 + Vl +4b : k = 0, ±1, ±2, ±3, ±4, ±5, ±6 L m G {4, 6, 12}. 



Step 2. Find D ^% defined by 
D = n £ Z + : n 



In this case we start with m = 4 to obtain the solution, afterwards m = 6 and 
finally m = 12. IfL> = 0, then the differential equation is not integrable because 
it falls in the case 4. Now, if Card(D) > 0, then for each n6fl with its respective 
to, search a rational function 

12^ x-c 

cer' 

and a polynomial S defined as 

s=n(x-c). 

cer' 

Step 3. Search for each n G D, with its respective m, a monic polynomial 
P n = P of degree n, such that its coefficients can be determined recursively by 

P_i = 0, P m = -P, 



P_i = -55 x p - ((to - i) d x S - SO) p - (to -%){%+ 1) S 2 rP i+1 , 

where i G {0, 1 . . . , to — 1, m}. If P does not exist, then the differential equation 
is not integrable because it falls in Case 4. Now. if P exists search w such that 

i=0 y ' 

then a solution of the differential equation is given by 
where u> is solution of the previous polynomial of degree to. 

Remark 1.1.14 ([5]). If the differential equation falls only in the case 1 of Kovacic's 
algorithm, then its differential Galois group is given by one of the following groups: 
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11 e when the algorithm provides two rational solutions. 

12 G' n l when the algorithm provides two algebraic solutions Ci , (,2 such that £™i Ca G 



C" S C(x) with n minimal. 

14 G m when the algorithm provides two non-algebraic solutions. 

15 <G a when the algorithm provides one rational solution and the second solution 

is not algebraic. 

16 B when the algorithm only provides one solution £ such that £ and its square 

are not rational functions. 

Kovacic's Algorithm in Maple 

In order to analyze second order linear differential equations with rational co- 
efficients, generally without parameters, a standard procedure is using Maple, 
and especially commands dsolve and kovacicsols. Whenever the command 
kovacicsols yields an output "[ ]", it means that the second order linear dif- 
ferential equation being considered is not integrable, and thus its Galois group is 
non- virtually solvable. 

In some cases, moreover, dsolve makes it possible to obtain the solutions in 
terms of special functions such as Airy functions, Bessel junctions and hypergeo- 
metric functions, among others (see [1]). 

There is a number of second order linear equations whose coefficients are 
not rational, and whose solutions Maple can find with the command dsolve but 
the presentation of the solutions is very complicated, furthermore the command 
kovacicsols does not work with such coefficients. These problems, in some cases, 
can be solved by our algebrization method (see section T2 .3. II and see also [5]). 

Beyond Kovacic's Algorithm 

According to the works of Michael Singer and Felix Ulmer [Ml [Ml [Ml [Ml QUO], 
we can have another perspective of the Kovacic's algorithm by means of the m-th 
symmetric power of a linear differential equation C, which is denoted as 

Theorem 1.1.15. Let C be a linear homogeneous differential equation of arbitrary 
order n. For any m > 1 there is another linear homogeneous differential equation, 
denoted by with the following property. If £1 , . . . , £ n are any solutions of C 



then any homogeneous polynomial in £1 j • • • , Cn of degree m is a solution of C^ m . 
In this way, Kovacic's algorithm can be stated as follows (see also [33j [34j 





[TOT]). 
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Algorithm. Consider C = d 2 ( — r(. 

Step 1. Check if C has an exponential solution £ = e-*" (with u € C(.x)). If so 
then C = ef u . 

Step 2. Check if C® 2 has an exponential solution £ = e-f 11 (with it G C(x)). If so, 
let v be a root of 

a 1 2 A 

v +uv+l-o x u+-u —u\=0. 

Then ( = e-K 

Step 3. Check if C® 4 , £® 6 , or C® 12 has an exponential solution Q = eJ u (with 
u € C(.t)). If so, then there is a polynomial of degree 4, 6 or 12 (respectively) 
such that if v is a solution of it then ( = e-> v . 

This algorithm can be generalized using the results of Michael Singer in [91] . 
The trick is to find the correct numbers (like 2, 4, 6, 12 of the Kovacic's algorithm). 

Theorem 1.1.16. Suppose a linear homogeneous differential equation of order n is 
integrable. Then it has a solution of the form 

C = eJ v 

where v is algebraic overC(x). The degree of v is bounded by I{n), which is defined 
inductively by 

1(0) = 1 

I(n) = max{ J(n), n\I(n — 1)} 



2rT / s 2rT 



J(n) = (VSn + l)' 

The following theorem is due to Singer and Ulmer (see |95jV 

Theorem 1.1.17 (Singer & Ulmer, [95j). Let C be a homogeneous n-th order linear 
differential equation over K with differential Galois group DG&l(L/K) C GL(n, C). 

1. If C has a Liouvillian solution whose logarithmic derivative is algebraic of de- 
gree m, then there is a DGal(L/ K) -semi-invariant of degree m inC\Y\, . . . , Y n ] 
which factors into linear forms. 

2. If there is a DGal(L/ K) -semi-invariant of degree m in C[Y±, . . . ,Y n ] which 
factors into linear forms, then C has a Liouvillian solution whose logarithmic 
derivative is algebraic of degree < m 

This algorithm is not considered implementable; the numbers I(n) are simply 
much too large. For example 1(2) — 384064, so this algorithm would require 
checking if C® m has a solution C with d x (/( e C(x) for m = 1,2,..., 384064. 
However, specific algorithms exist for order three and higher, see [40l [95) [96l [99l 
fTOOl 146], 



1.1. Picard-Vcssiot theory 



21 



1.1.3 Eigenrings 

We consider two different formalisms for Eigenrings, the matrix and operators for- 
malism. We start with the Matrix formalism of Eigenrings following M. Barkatou 
in [TJ, but restricting again to 2 x 2 matrices. 

Let K be a differential field and let A be a matrix in GL(2, K) such that, 

d x X = -AX. (1.8) 

Consider a matrix equation (jTTHJ) and let P 6 GL(2, K). The substitution X = PY 
leads to the matrix equation 

d x Y = -BY, B = P~ l {d x P + AP). (1.9) 

Definition 1.1.18. The matrices A and B are equivalents over K, denoted by A ~ 
P, when there exists a matrix P £ GL(2, if) satisfying equation (|1.9p . The systems 
(jl.8p and (jl.9[) are equivalent, denoted by [A] ~ [P], when A and P are equivalent. 

By equation (|1.9[) . we have PP = c^P + AP. In general, assuming PP = 
PA, where P is a 2 x 2 matrix, i.e, P is not necessarily in GL(2,K), we obtain 
PA = d x P + AP, which leads us to the following definition. 

Definition 1.1.19. The Eigcnring of the system [A], denoted by £(A), is the set of 
2x2 matrices P in K satisfying 

d x P = PA- AP. (1.10) 

Equation fTTTUj) can be viewed as a system of 4 first-order linear differential 
equations over K. Thus, £(A) is a C— vector space of finite dimension < 4. Owing 
to the product of two elements of £(A) is also an element of £{A) and the identity 
matrix I2 belongs to £(A), we have that £ (A) is an algebra over C, i.e., £(A) is a 
C— algebra. As a consequence, we have the following results that can be found in 
0- 

Proposition 1.1.20 (0). Any element P of£(A) has: 

• a minimal polynomial with coefficients in C and 

• all its eigenvalues are constant. 

Proposition 1.1.21 (|7J). If two systems [A] and [B] are equivalent, their eigenrings 
£{A) and£{B) are isomorphic as C— algebras. In particular, one has dimc£(A) = 
dim c £(P). 

Definition 1.1.22. The system [A] is called reducible when A ~ P, being P given 

by 

B =t b 

\021 022 

When [A] is reducible and 621 = 0, the system [A] is called decomposable or 
completely reducible. The system [A] is called irreducible or indecomposable when 
[A] is not reducible. 
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Assume that the eigenring £(A) is known. 

Theorem 1.1.23 ([7]). If £{A) is not a division ring then [A] is reducible and the 
reduction can be carried out by a matrix P G GL(2,K) that can be computed 
explicitly. 

The condition £(A) is not a division ring implies dime £ (A) > 1. Indeed, if 
P G £{A) \ {0} is not invertiblc, then the family {/, P} is linearly independent 
(over C) and hence dime £(A) > 1. In our case the converse is true, due to the 
field of constants C is algebraically closed. Indeed, suppose that dime £ (A) > 1 
then there exists P G £(A) such that the family {I,P} be linearly independent. 
Since C is algebraically closed, there exists A G C such that det(P — XI) = 0. Hence 
£(A) contains an element, namely P — XI, which is non-zero and non invertible. 

The computation of eigenrings of the system [A] is implemented in ISOLDE 
(Integration of Systems of Ordinary Linear Differential Equations). The function 
is eigenring, the calling sequence is eigenring (A, x) being the parameters: A - 
a square rational function matrix with coefficients in an algebraic extension of the 
rational numbers and x - the independent variable (a name) . ISOLDE was written 
in Maple V and it is available at http://isolde.sourceforge.net/. 

In operators formalism, we restrict ourselves to second order differential op- 
erators and we follow the works of Singer, Barkatou and Van Hoeij (see [93l [3 HH 
SH [45]). A differential equation C := d x y + ad x y + by = with a, b G K corre- 
sponds to a differential operator / = d\ + ad x + b acting on y. The differential 
operator / is an element of the non-commutative ring i^ffij,]. 

The factorization of operators is very important to solve differential equa- 
tions, that is, a factorization / = ££K where £, £H G A'[9 X ] is useful for computing 
solutions of / because solutions of the right-hand factor £H are solutions of / as 
well. 

Definition 1.1.24. Let £ be a second order differential operator, i.c C := £(y) = 0. 
Denote V(£) as the solution space of C. The Eigenring of £, denoted by £(£), is 
the set of all operators *K for which 9t(V(£)) is a subset of V(£), that is £5H = 6£, 
where 6 is also an operator. 

As consequence of the previous definition, 9\ is an endomorphism of the so- 
lution space V(£). This means that we can think of 91 as a linear map V — > V and 
choosing one local basis of V we obtain, by linear algebra, that £H has a matrix 
Mjh. The characteristic polynomial of this map can be computed with the classi- 
cal methods of linear algebra. For endomorphisms £H, the product of £ and D\ is 
divisible on the right by £. This means that if £(y) = 0, then £(9t(y)) = 0, so 
that m map V -> V. 
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For the general case of operators 



£ = J2 b kd*, m=J2 a i d ii ai,heK, n>m, C := £(y) = 

fc=0 i=0 

and G = DGal(L/K), we can see £H as a G-map. Now, denoting by P the char- 
acteristic polynomial of £H, assume that there exists polynomials Pi, P2 with 
gcd(Pi,P2) = 1 such that P = P1P2. By Cayley-Hamilton theorem we have that 
P(Mm) = and by kernel theorem we have that V = ker(Pi(M<K))©ker(P 2 (MfR)), 
in where ker(Pi(Aifn)) and ker(p2(M(n)) are invariants under G. 



Let A denote an eigenvalue of 9\; there exists a non-trivial eigenspace V\ C 
V, which means that £ and 9\ — A has common solutions and therefore right — 
gcd(£, £H — A) is non-trivial, this means that right — gcd(£, 9\ — A) it is a factor of 
£. 

Returning to the second order operators, we establish the relationship be- 
tween the Eigcnring of the system [A] and the Eigcnring of the operator £. We 
start recalling that C given by d^y + ad x y + by — 0, a, b G K, can be written as 
the system 

' J ) = (\ 1 ) ( B 
dxV J \ -b -a J \ d x y 



Or 

and the system of linear differential equations 



z J \ c d J \ z 
by means of an elimination process, is equivalent to the second-order equation 

dly- (a + d+^-jd x y- (d x a + be- ad- a^j y = 0, b^O. (1.11) 

In this way, we can go from operators to systems and reciprocally computing 
the Eigcnrings in both formalism. In particular, we emphasize in the operator 
£ = d x + pd x + q, which is equivalent to the system [A], where A is given by 

As immediate consequence we have the following lemmas. 
Lemma 1.1.25. Consider £, A and P as follows: 

£ = dl+pd x + q, A =(° q > P= (c d) ' a,b,c,d,p,qe K. 

The following statements holds 
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1. IfPe 8(A), then = a + bd x € 8(2). 

2. Ify\ = a + bd x e 8(2), then P e £ (A), w/iere P is given by 



P = 



a b 
d x a — bq a + 9^6 — bp 



3. 1 < dim c £(£) < 4. 
4- PGGL(2,Jr)^^f -f + p^ 
Proof. Suppose that Vi and Va are solution spaces of C and [A] respectively. 

1. If P G 8(A), then for a solution Y = (^j ofd x Y = -AY, if we set Z = PY 
then d x Z = —AZ. Setting 

we have 

d x (i = Ca, 

Ci = ay + 69 x y. 

Owing to £(Ci) = for all y G Vc, a + foc^ maps to Vjj and so 9t := 
a + 6^ G 8(2). 

2. U m = a + bd x e 8(2), then let C := W(y) = ay + bd x y. Then 2(VH(y)) = 
2(C) = and 



Thus, there exists P satisfying 



t \ - p( y 



Jd x Cj \<>.<l) 1 

which lead us to the expression 
C = ay + bd x y, 

d x ( = (d x a - bq)y + (a + bd x - bp)d x y. 
Statements 3 and 4 can be obtained immediately. □ 



Remark 1.1.26. Using this P, one can obtain the characteristic polynomial of P 
which is the characteristic polynomial of 9t and so obtain eigenvalues of 91 needed. 
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Remark 1.1.27. Let £ be the differential operator d x + b, where b £ K, C := 
£(y) = 0. The dimension of the eigenring of £ is related with: 

• the number of solutions over K of the differential equation C and its second 
symmetric power and 

• the type of differential Galois group (see [7] l93 l l43 j l44l [45] ) . 

The previous remark is detailed in the following lemma. 

Lemma 1.1.28. Assume £ = d\ + b, where b € K , L := £(y) = 0. The following 
statements holds. 

1. //dimc£(£) = 1, then either differential Galois group is irreducible (Hoo or 
primitive), or indecomposable (G CM, G ^ {e, G m , G a , G^, G^}). 

2. J/dimc£(£) = 2, then either, the differential Galois group is the additive 
group or is contained in the multiplicative group, but never is the identity 
group. Thus, we can have two solutions but not over the differential field K . 

3. J/dimc£(£) = 4, then the differential Galois group is the identity group. In 
this case we have 2 independent solutions Ci and (2 in which Cf an d C1C2 
are elements of the differential field K , i.e. the solutions of belongs to 
the differential field K . 

Proof. Assume that G = T)G&\(L/K), then G commutes with each element P G 
£(£), i.e., GP = PG. Computing the dimension of the set of matrices satisfying 
GP = PG we obtain dim c £(£) for each G C SL(2, C). □ 

The eigenring for a differential operator £ has been implemented in Maple. 
The function is eigenring, the calling sequences are eigenring(L, domain) and 
endomorphism-charpoly (L, R, domain) , where £ is a differential operator, *H is 
the differential operator in the output of eigenring. The argument domain describes 
the differential algebra. If this argument is the list [Dx,x] then the differential 
operators are notated with the symbols Dx and x, where Dx is the operator d x . 

Example. Consider K = C(x), £ = d\ — ■% and [A] where A is given by 

A = {-6x- 2 0) ■ 
The Eigenring of [A] and the Eigenring of £ are given by 

£(A)=Vect((j i)\ 6 ~-i o)'(-to 4 3^) '(-4a- 8 -2a- 5 ))' 

£(£) =Vect( l,xd x -l,x 6 d x -3x 5 ,^r + ) . 

\ nr. r ^ I 
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1.1.4 Riemann's Equation 

The Riemann's equation is an important differential equation which has been 
studied for a long time, since Gauss, Riemann, Schwartz, etc., see for example 
[521 ITS] . We are interested in the relationship with the Picard-Vessiot theory. 
Thus, we follow the works of Kimura [54], Martinet & Ramis [64] and Duval & 
Loday-Richaud [55] , 

Definition 1.1.29. The Riemann's equation is an homogeneous ordinary linear dif- 
ferential equation of the second order over the Riemann's sphere with at most 
three singularities which are of the regular type. Assuming a, b and c as regular 
singularities, the Riemann's equation may be written in the form 

„o ( \ — p — d 1 — a — a' 1 — r — t j \ „ 

d 2 x V + fLJL + r— + d x y 1.12 

\ x — a x — b x — c ) 

+ ( pp'{a-b)(a-c) + ao->(b-a)(b-c) + tt> (c - a)(c - b) \ = Q 
\ (x — a) 2 (x — b)(x — c) (x — b) 2 {x — a){x — c) (x — c) 2 (x — a)(x — b) J 

where {p,p'), (a, a 1 ) and (r, r') are the exponents at the singular points a, b, c 
respectively and must satisfy the Fuchs relation p + p' + a + a' + t + t' = 1. 
The quantities p' — p, a' — a and r' — r are called the exponent differences of the 
Riemann's equation (|1.12| at a, b and c respectively and are denoted by A, p and 
v as follows: 

A = p — p, p, = a' — a, v = t — t. 

The complete set of solutions of the Riemann's equation (|1.12p is denoted by the 
symbol 

a b c 
y = P { p a t x 

p' G 1 T' 

and is called Riemann s P -function. 

Now, we will briefly describe here the Kimura's theorem that gives neces- 
sary and sufficient conditions for the integrability of the Riemann's differential 
equation. 

Theorem 1.1.30 (Kimura, [54]). The Riemann's differential equation (|1.12p is in- 

tegrable if and only if, either 

(i) At least one of the four numbers X + Ji + v, —A + ju + v, X — Jl + v, X + p — v 
is an odd integer, or 

(ii) The numbers X or —A, p or —p and v or —v belong (in an arbitrary order) 
to some of the following fifteen families 
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1 

1 


1 /O i 7 

1/2 + t 


1/Z + 771 


arbitrary complex number 




2 




l/Q | „-„ 
1/0 + 777 


I/O 1 „ 

1/3 + q 




o 
o 


9 /Q _i_ 7 

2/o-r l 


1/ O -|- 777 


1/6 + q 


— 

/ + 777 + 9 even 


4 


1/2 + I 


1/3 + 772 


1/4 + q 




5 


2/3 + 1 


1/4 + 777 


1/4 + g 


I + ?77 + 9 even 


6 


1/2 + 1 


1/3 + 777 


1/5 + g 




7 


2/5 + 1 


1/3 + 777 


1/3 + g 


I + m + 9 even 


8 


2/3 + / 


1/5 + 777 


1/5 + 9 


/ + ?77 + 9 even 


9 


1/2 + / 


2/5 + 777 


1/5 + g 


/ + m + 9 even 


10 


3/5 + / 


1/3 + 777 


1/5 + g 


/ + 777 + 9 even 


11 


2/5 + / 


2/5 + 777 


2/5 + q 


/ + m + 9 even 


12 


2/3 + / 


1/3 + 777 


1/5 + 9 


/ + ?77 + 9 even 


13 


4/5 + / 


1/5 + 777 


1/5 + 9 


/ + m + 9 even 


14 


1/2 + / 


2/5 + 777 


1/3 + 9 


/ + m + 9 even 


15 


3/5 + / 


2/5 + 777 


1/3 + 9 


/ + m + 9 even 



/fere n, m, 9 are integers. 

Using the Moebius transformation [41) . also known as nomographic substi- 
tution, in the Riemann's equation (jl,12[) . we can map x = a,b,c to x' = a' , b' , d , 
respectively: 

, px + q 
x = . 

rx + s 

In particular, we can place the singularities at x — 0, 1, 00 to obtain the following 
Riemann's equation: 



1 - p - p' l-a-a' 



PP 

„2 



x - 1 
aa' 



{x-iy 2 



d x y 



pp era 



x(x — 1) 



(1.13) 



y = o, 



where the set of solutions is 



y = P 




Sometimes it is very useful to map x = 0, 1, 00 to x' = —1, 1, 00 in the Riemann's 
equation (j!.13[) . for example, setting p = 0, we can state the substitution: 




-1 1 00 
a a 2t y/x 
a' a 1 2t' 



We can transforms equation (| 1 . 13[) to the Gauss Hypergeometric equation as 
follows: 
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1 oo 

x p (x- l) a P { k x 

1— 7 7 — k — 

where k = p + <j + t, f3 = p + o~ + t' and 7 = 1 + p — p' . Then 

( 1 00 

y = P< k x 

I 1 — 7 7 — k — 

is the set of solutions of the Gauss Hypergeometric differential equation^ 

^ + 1(1- ,) ^ - - °< 

where the Fuchs relation is trivially satisfied and the exponent differences are given 

A = 1 — 7, fx = 1 — 7 — 0, V = — K. 

We remark that the Galoisian structure of the Riemann's equation do not change 
with the Moebius transformation. 

The confluent Hypergeometric equation is a degenerate form of the Hyper- 
geometric differential equation where two of the three regular singularities merge 
into an irregular singularity. For example, making "1 tend to 00" in a suitable way, 
the Hypergeometric equation (11.141) has two classical forms: 



Kummer's form 



• Whittaker's form 



%y + ?—Zd aV -±y = (1.15) 



where the parameters of the two equations are linked by k = | — a and fj, = § — 5 • 
Furthermore, using the expression (|1.3p . we can see that the Whittaker's equation 
is the reduced form of the Kummer's equation. The Galoisian structure of these 
equations has been deeply studied in [MJ [35] . 

Theorem 1.1.31 (Martinet & Ramis, [MJ). The Whittaker's differential equation 
(11.16[) is integrable if and only if either, k + p E | + N, or k — p E 5 + N, or 
-K + p, S \ + N, or -K - fj, e 5 + N. 



1 In general, for the Hypergeometric differential equation, is used a instead of k, but we want 
to avoid further confusions. 
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The Bessel's equation is a particular case of the confluent Hypergeometric 
equation and is given by 

%y+-d x y+?—^-v = 0. (1.17) 
x x z 

Under a suitable transformation, the reduced form of the Bessel's equation is a 
particular case of the Whittaker's equation. Thus, we can obtain the following well 
known result, see [55] p. 417] and see also [5rJI IrJT] . 

Corollary 1.1.32. The Bessel's differential equation (|1.17[) is integrable if and only 

ifn e 5 +Z. 

We point out that the integrability of Bessel's equation for half integer of 
the parameter was known by Daniel Bernoulli [104] . By double confluence of the 
Hypergeometric equation (|1 . 14[) , that is making "0 and 1 tend to oo" in a suitable 
way, one gets the parabolic cylinder equation (also known as Weber's equation): 

bty=(jx 2 - l --n^jy, (1.18) 

which is integrable if and only if n G Z, see [56] [33]. Setting n — ~ \ and 

making the change x i— > ^J~^(ax + 6), one can gets the Rehm's form of the Weber's 
equation: 

dly = (ax 2 + 2bx + c) y, a ^ 0, (1.19) 

so that - — - is an odd integer. 

The Hypergeometric equation, including confluences, is a particular case of 
the differential equation 

d ly + ^ d *y + ^y> AgC ' L = a + a lX , Q = b + b lX + b 2 x 2 . (1.20) 

We recall that the classical orthogonal polynomials and Bessel polynomials are 
solutions of equation (TOO"]) , see [2"31 ICTl 170] : 

• Hermite, denoted by H n , 

• Chebyshev of first kind, denoted by T n , 

• Chebyshev of second kind, denoted by U n , 

• Legendre, denoted by P n , 

• Laguerre, denoted by L n , 

• associated Laguerre, denoted by 

• Gegenbauer, denoted by ci 
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• Jacobi polynomials, denoted by Vn' v and 

• Bcssel polynomials, denoted by B n . 

In the following table we give Q, L and A corresponding to equation (|1 .20|) 
for classical orthogonal polynomials and Bessel polynomials. 



Polynomial 


Q 


L 


A 


H n 


1 


-2x 


2n 


T n 


l-x 2 


—x 


n 2 


u n 


1-x 2 


—3x 


n(n + 2) 


Pn 


l-x 2 


-2x 


n(n + 1) 


L n 


X 


l-x 


n 


T (rn) 


X 


m + I — x 


n 


r<( rn ) 


l-x 2 


-(2m + l)x 


n(n + 2m) 


In 


l-x 2 


v — m — (m + v + T)x 


n(n + 1 + m + v) 


B n 


x 2 


2(3 + 1) 


-n(n + 1) 



The associated Legendre polynomials, denoted by Pi m \ does not appear in 
the previous table. They arc solutions of the differential equation 



2x /n(fi+l)-^\ 

This equation can be transformed into a Ricmann's differential equation through 
the change x i— > 1 _ 1 a;2 . Thus, the complete set of solutions of equation (|1.2ip is 
given by 




the exponent differences are A = h, Jl = | and v — 0. By Kimura's theorem this 
equation is integrable. 



Finally, we remark that integrability conditions and solutions of differential 
equations with solutions orthogonal polynomials, including Bessel polynomials, 
can be obtained applying Kovacic's algorithm. In the same way, we can apply 
Kovacic's algorithm to obtain the same results given by Kimura [S3] and Mar- 
tinet & Ramis [64] . Also we recall that Duval & Loday-Richaud applied Kovacic's 
algorithm to some families of special functions [33] . 
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1.2 Supersymmetric Quantum Mechanics 

In this section we establish the basic information on Supersymmetric Quantum 
Mechanics. We only consider the case of non-relativistic quantum mechanics. 

1.2.1 The Schrodinger Equation 

In classical mechanics the Hamiltonian corresponding to the energy (kinetic plus 
potential) is given by 

II~P|| 2 -> -> -> 

H = — \-U(x), p = (pi,...,p n ), X =(x 1 ,...,X n ), 

2m 

while in quantum mechanics the momentum ~p is given by p = — ifiV, the Hamil- 
tonian operator is the Schrodinger (non-relativistic. stationary) operator which is 
given by 

and the Schrodinger equation is = , where x is the coordinate, the eigen- 
function 'J is the wave function, the eigenvalue E is the energy level, U(x) is the 
potential or potential energy and the solutions of the Schrodinger equation are the 
states of the particle. Furthermore, it is known that = H, i.e., the Schrodinger 
operator is a self-adjoint operator in a suitable complex and separable Hilbert space. 
Thus, H has a purely real spectrum spec(_ff) and its spectrum spec(_ff ) is the dis- 
joint union of the point spectrum sj>ec p (H) and the continuous spectrum spec c (iJ), 
i.e., spec(-ff) = spec p (iJ) Uspec c (-ff) with spec p (H) n spec c (iJ) = 0. See for exam- 
ple 0EE1G1]. 

Along this thesis we only consider the one-dimensional Schrodinger equation 
written as follows: 

H^ = E^, H = -d 2 z + V(z), (1.22) 

where z — x (cartesian coordinate) or z = r (radial coordinate) and h = 2m = 1. 
We denote by ^ n the wave function for E = E n . The potentials should satisfy 
some conditions depending of the physic situation such as barrier, scattering, etc., 

see [mi Eg eqi Eg ng. 

Definition 1.2.1 (Bound States). The solution ^„ is called a bound state when E 
belongs to the point spectrum of H and its norm is finite, i.e., 

E„ G spec p (iJ), J \^ n (x)\ 2 dx < oo, n G Z+. (1-23) 

An interesting property of bound states is given by the Sturm's theorem, see 

El Eg. 
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Theorem 1.2.2 (Sturm's Theorem). If • • • > ^n, ■ ■ ■ are the wave functions 

of the bound states with energies E < E\ < ■ ■ ■ < E n < . . ., then has n 
nodes (zeros). Furthermore, between two consecutive nodes of^ n , there is a node 
of and moreover ^ n + r has at least one zero for all r > 1. 

Definition 1.2.3 (Ground State and Excited States). Assume "fo, ^lj • • • s ^n, ■ ■ ■ 

as in the Sturm's theorem. The state ^Sq, which is state with minimum energy is 
called the ground state and the states . . . , \£„, . . . are called the excited states. 

Definition 1.2.4 (Scattering States). The solution iff corresponding to the level 
energy E is called a scattering state when E belongs to the continuous spectrum 
of H and its norm is infinite. 

The wave function belonging to the continuous spectrum have two typical 
boundary conditions: the first ones, barrier potentials and the second one periodic 
boundary conditions. The transmission and reflection coefficients are related with 
the barrier potentials, |36j . 

When the particle moves in one dimension, we use the classical one dimen- 
sional Schrodingcr equation with cartesian coordinate x. 

Example (The Harmonic Oscillator). We consider the Hamiltonian operator H 
given by (|1.22|) . Normalizing the angular velocity (lu = 1), the one- dimensional 
harmonic oscillator potential is V = \x 2 . The creator (raising) and annihilator 
(lowering) operators, denoted respectively by a) and a, given by 

a f = -d x + a = d x + \<x 

lead us to the relations 

H = ata+-, [0,0^=1, [at,B\ = -at, [a,H]=a. (1.24) 

We want to solve H^ff = E^ff, in particular we are interested in the case 
Hiff n = E n ^ n , where is the ground state and ^f\,...,^f n are the excited 
states. Considering H^ffo = Eo'ffo and assuming that a'a^ffo = 0, we obtain the 
energy Eq and the ground state \fo : 

Using the relation *ff n = (a') n ^fo, we obtain the rest of wave functions 
^ = e -i x \ *! = at* = xe~i x \ ^ 2 = a^ 1 = (x 2 -l)e-i x \ 



* 3 = {x 3 - 3.t)* , *4 = (x 4 -6.t 2 + 3)*o,--- *n = a^ n -i = H n *o, 
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where H n are the Hermite's polynomials presented in the previous section. Now, 
to obtain the complete energy spectrum we use the relations (|1.24p , thus 

1 3 5 2n + 1 

— 7j' Ei — I + E — -, E 2 = 1 + -Ei = -, . . . , E n — — - — . 

Another way to obtain E n is given by the formula 

E n = ^- d J^. 
4 * n 

Also we can consider the particle moving in three dimensions, this means 
that ~p = {p x ,Py,p z ), where p x = -td x , p y = -id y , p z = -id z and p = The 
angular momentum operator is given by L = {L x , L y , L z ) where L x = yp z — zp y , 
L y = zp x — xp z , L z = xp y — yp x and L = \\L\\. The square of the angular 

momentum operator || L || 2 = L 2 = L\ + L y + L 2 Z commutes with all components 
of the angular momentum operator. 

In spherical coordinates x ~ r sin v cos ip, y = rsin^sin^, z = rco&v, L 2 is 
given by 

L 2 = -A v . v , A»», = -J—d u (smvd,,) + \ d 2 
sin v sin v v 

where we denote by that part of the Laplacian acting on the variables v and 
(p only. The kinetic energy given by p 2 = d 2 + d 2 + d z reads in polar coordinates 
as 

p 2 = P 2 r + ^L 2 , pl = ~dr^d r )=-(d 2 r + 2 -d r 

Now, for central potentials, where the potential U^r) is spherically symmetric, 
i.e., U( r ) = U(r), we can reduce the Schrodinger equation to an one dimensional 
problem, the so-called radial equation. 

We start writing the eigenfunctions and eigenvalues of the operator L 2 : 

L 2 Yu n {v, i P )=l{t+l)Y l , m {v, if), 

the eigenfunctions Y(y, ip) are the spherics harmonics which are related with the 
associated Legendrc Polynomials 

Yt, m (v,<p) = PZ n {cosvy m *. 

Assuming $ as eigenfunctions of p 2 + U(r) satisfying $ = J?^(r)l^ jm (f, cp), 
i.e., the partial wave function decomposition sec ([36, 60, 88 J, we have 



Pr + ^2 L + u ( r ) - E ) Rt(r)Yt, m (v, ip) = 0, 
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so that we obtain 

(rf + U(r) - E) *(r) + ^^fY^, *) = 

and owing to L 2 Y^ m (z/, 95) = 1(1 + l)Yi tJn [y, if) we have the radial equation 

pi + £ JI±}1 + u( r )\ R e {r) = ER e (r). 

Applying the expression (|1.3p , the radial equation can be reduced to the Schrodinger 
equation (| 1 . 22[) as follows: 

H^ = E^ 1 H = -$ + V(r), V(r) = i JA±ll + jj( r ), V = rR t (r). 

The equation for the angular part is always solved through spherics harmon- 
ics, while for the radial part, the analysis depends on the spherically symmetric 
potential U(r). One example of the radial equation is the Coulomb potential. The 
complete set of physical and mathematical conditions for the potentials, spec- 
trum and wave functions in one or three dimensions can be found in any book of 
quantum mechanics, including bound states and scattering cases, see for example 

MMM- 



1.2.2 Darboux Transformation 

The following theorem is the most general case for Darboux transformation in the 
case of second order linear differential equations, which is taken faithfully from 
[28]. 

Theorem 1.2.5 (Darboux, [28]). Suppose that we know how to integrate, for any 
value of the constant m, the following equation 

d 2 x y + Pd x y+{Q-mR)y = Q. (1.25) 

If 6 7^ is an integral of the equation 

d 2 J + PdJ + Q6 = 0, 

then the function 

d x y - ^-y 

u= — 8 \ (1.26) 



IR 

will be an integral of the equation 

d\u + Pd x u + [e^R.d x (^=j - eVRd 2 x \^=) 

for m =/= 0. 
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Darboux in [281 [29] presented the particular case for R = 1 and P = 0, which 
today is known as Darboux transformation, but really is a corollary of the general 
Darboux transformation given in theorem 11.2.51 

Corollary 1.2.6 (Darboux. [281 129j ). Suppose that we know how to integrate 

d 2 x y = (f(x)+m)y (1.28) 
for any value of m. If 6 satisfies the equation d 2 9 = (f(x) + m{)6, the function 

u = a x y —y 



will be an integral of the equation 

dlu = (e&i | 

for m ^ mi. Furthermore, 



1\ / d x 8\ ( d x 



2 



0d 2 x - rni = fix) -2d x \^-j=2\^-j - f{x) - 2mi. 

Remark 1.2.7. In practice, we need two values of m to apply the Darboux's results. 

Example. Consider the equation d x y = my . Employing the solution 9 = x, we 
shall get 

d ly = f~-r + TO ) V- 

Applying the same method to the latter equation, but taking now 6 = x 2 , we shall 
have 

d *y = I ~ + m J y 

and so on. The cases mi = and mi = — 1 also can be found as exercises in the 
Ince's book (HI p. 132]. 

We can see that equation (|1.28[) coincides with the Schrodinger equation 
(jl.22[) . Thus, we can apply the Darboux transformation in where m = — E and 
mi = -E Q . 

The following definition corresponds with Delsarte's transformation opera- 
tors (isomorphisms of transmutations), which today are called intertwiner opera- 
tors, see [30] . 

Definition 1.2.8. Two operators £o and £i are said to be intertwined by an oper- 
ator T if 

£iT = T£ - (1-30) 
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We can relate the intertwiner operators with the Darboux transformation of 
equation p.22[) . where £1 and £ are Schrodinger operators and T can be either 

+ x ^ *0 ' 

Crum, inspired by the works of Liouville [61 , 62J obtained one kind of iterative 
generalization of Darboux's result giving emphasis in the Sturm-Liouvillc systems, 
i.e., he proved that the Sturm-Liouville conditions are preserved under Darboux 
transformations, see [27] . The Crum's result is presented in the following theorem, 
defining the Wronskian determinant W of k functions fx, fi, ■ ■ ■ , fk by 

W(fx,...,f k )=detA, A,, iK i.j 1.2 k. 

Theorem 1.2.9 (Crum, [27]). Let ^x,^2, ■ ■ - ^n be solutions of the Schrodinger 
equation (|1.22p for fixed, arbitrary energy levels E = Ex , E<i , ■ ■ ■ , E n , respectively. 
Then, we obtain the Schrodinger equation 

ffW*[n] =EV[n], E^Ei,l<i<n, H [n] = -8% + V[n], 

where 

m = V[n] =V- 2dl In W(*x, *«)• 

Darboux transformation coincides with Crum's result in the case n = 1 and 
the iterations of Darboux transformation coincides with Crum iteration, see [71j . 
Both formalisms allow us to obtain new families of Schrodinger equations preserv- 
ing the spectrum and the Sturm-Liouville conditions, see |66[ 171]. Furthermore, 
there are extensions of Crum's iteration connecting the Sturm-Liouville theory 
with orthogonal polynomial theory |58j . 

Schrodinger in [89] factorized the Hypergeometric equation (|1 . 14|) . He started 
making the change of variable 2x — 1 = cos 8, after, using the expression (|1.3p . he 
reduced the Hypergeometric equation to obtain conditions of factorization. In this 
way, setting k/3 = E, we obtain families of Schrodinger equations Q1.22p . This re- 
sult was used by Natanzon in [69j to obtain the well known Natanzon's potentials, 
i.e, potentials which can be obtained by transformations of the Hypergeometric 
equation and its confluences, see [24\ 125] . In particular, the Ginocchio potentials 
are obtained through the Gegenbauer polynomials. 

Witten in [1071 §6] presented some models in where dynamical breaking 
of supersymmetry is plausible. The first model is a model in potential theory- 
supersymmetric quantum mechanics, although is not a model in the field theory. 

Definition 1.2.10. A super symmetric quantum mechanical system is one in which 
there are operators Qi, that commute with the hamiltonian TL, 



[Qi,H]=0, i = l,...,n 



(1.31) 
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and satisfy the algebra 

{Q 4 , Qj} = StjH, {Q u Qj} = Q l Q J + Q J Q l . (1.32) 

The simplest example of a supersymmetric quantum mechanical system cor- 
responds to the case n = 2, which in physical sense involves a spin one half particle 
moving on the line. This case is the main object of this thesis. The wave function 
of 7i$ = _E$ is therefore a two-component Pauli spinor, 



The supercharges Qi arc defined as 
aip±a 2 W(x) 

Q± = 5 ' Q+ = QuQ-=Q2, P = -io x , (1.33) 

where the superpotential W is an arbitrary function of a; and o~i are the usual Pauli 
spin matrices 

l\ {0 -i\ (\ 



: \\ 0j> <72= [i 0)> as= \0 -1 
Using the expressions (|1.31|) . (|1.32p and (|1 .33|) we obtain H. 

The supersymmetric partner Hamiltonians H± are given by 

The potentials V± are called supersymmetric partner potentials and arc linked 
with the superpotential through a Riccati equation. So that equation (|1.34|) 
can be written as 

'H+ 
PL- 



n = 



which lead us to the Schrodingcr equations = E'$> + and H-ty- = Efy-, 

and for instance, to solve 7i$ = EQ is equivalent to solve simultaneously H + ^ + = 
and = 

We analyze equation Qi& = 0, which must be satisfied by a supersymmetric 
state. This is due to Q\ = Q\ = ^TL and for instance Q2 = —icrsQi, which implies 
that <5i$ = 0, or aip<& = —a 2 W^>. Now, multiplying by o\ and using the facts 
that p = —id x , o~\&2 — io~3, this equation becomes 

d x § = W(x)o- 3 $(x), (1.35) 
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and the solution is 

In agreement with [107] , one important generalization of the previous model 
could be its extension to four dimensions, i.e., the case n = 4. 

V.B. Matveev and M. Salle in [66] interpret the Darboux Theorem as Dar- 
boux covariance of a Sturm-Liouville problem and proved the following result, sec 
also [H §5-6]. 



Theorem 1.2.11 (Matveev & Salle, [66j)- The case n = 2 in Supersymetric Quan- 
tum Mechanics is equivalent to a single Darboux transformation. 

According to Natanzon |69| . a solvable potential, also known as exactly solv- 
able potentials, is a potential in which the Schrodinger equation can be reduced to 
hypergeometric or confluent hypergeometric form. The following are examples of 
solvable potentials. 



Potential Name 

V(x) = I ^' X j E^' z\ Infinite square well 

I co, x & [U, Li\ 

V(r) = \ ' r , ' A Radial infinite square well (1-37) 

I oo, r $ [0,L] 

V(r) = —^hf, Hulthen 

V(r) = - eK t!'_ 1 + u^l i)3 1 Generalized Hulthen 

V{x) = \^x 2 +g aJ f^, g a >0 ^ 

We remark that our definition of intcgrability, definition 11.1.81 is different of 
the concept of solvability given by Natanzon. In the next chapter we come back 
on this problem. 

According to Dutt et al. |32j . a conditionally solvable potential is a potential 
in which the entire bound state spectrum can be analytically obtained, where the 
parameters in the potential satisfies a specific relation. The following potentials 
are two examples of conditionally solvable potentials 

v( , A_ _ __B 3 

[X> ~ 1 + e-^ Vl + e- 2 - 4(1 + e- 2 *) 2 

A Be' x 3 

{X) ~ 1 + e-^ ~ V e~ 2 * + 1 ~ 4(1 + e~ 2 *) 2 ' ( ° 8) 
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As a generalization of the method to solve the harmonic oscillator [37l [31] , 
the ladder (raising and lowering) operators are defined as 

A+ = A = d x -^, 

which are very closed with the supercharges Q± in the Witten's formalism. Thus, 
AV o = 0, A + A = H_, AA+ = H + = -dl + V + (x), where 

v + M-v. W -%(^)--v_ W + 

The supcrsymmctric partner potentials V+ and V— have the same energy levels, 
except for Eq ' = 0. In terms of the superpotential W(x), the operators A and 
A+ are given by 

A + = -d x + W(x), A = d x + W(x). 

In the same way, the supersymmetric partner potentials V± (x) and the su- 
perpotential W(x) satisfies: 

V+{X)+ 2 V - {X) =W*{ X ), [A,A+) = 2d x W{ X ). 

Let ^>n ' and \£„ ^ denote the eigenfunctions of the supersymmetric Hamilto- 
nians if_ and i7 + respectively, with eigenvalues i?i ' and E n + \ The integer 
n = 0, 1, 2, ... , denotes the number of nodes in the wave function. 



Theorem 1.2.12 (Dutt et al., [31]). If^h ' is any eigenfunction of H- with eigen- 
value En 
thermore 



value E„ ' , then A^ n ■* is an eigenfunction of H + with the same eigenvalue. Fur- 



£?(+) = E { ~ ] * (+ ) = A VI/ ( ~ ) 

Considering i?o = and n > in theorem 11.2.121 we can see that the su- 
persymmetric partner potentials V+ and K_ have the same spectrum. The ground 
state energy Eq = of V- has no corresponding level for V+. Furthermore, if the 
eigenfunction $L j of f/_ is normalized, then also the eigenfunction ^/ n + ^ of H+ 
is normalized. 

The operator A converts an eigenfunction of H_ into an eigenfunction of H+ 
with the same energy, whilst the operator A + converts an eigenfunction of H + 
into an eigenfunction of U_ with the same energy. Furthermore, the operator A 
destroys a node (^„+x has n+1 nodes, whilst \I>„ has n nodes) and the operator 
A + creates a node. In summary, the annihilator and creator operators (A and A + 
respectively) connect states of the same energy for two different supersymmetric 
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partner potentials. 

Gendenshtei'n, in his remarkable paper |37j . introduced the concept of shape 
invariance, which is a property or condition of some classes of potentials with 
respect to their parameters. Assuming a as a family of parameters, the shape 
invariancc condition can be seen such as follows: 

V n+1 (x; a n ) = V n (x; a n +i) + R(a„), V- = V , V + = Vi, 

where R is a remainder, which does not depends on x. In this way we say that 
V = V— = Vo is a shape invariant potential. 

The potentials, corresponding to a Schrodinger operator, satisfying the shape 
invariance property allow us to obtain a fully algebraic scheme for the spectrum 
and wave functions. This is illustrated in the following theorem obtained by Gen- 
denshtei'n. 

Theorem 1.2.13 (Gendenshtei'n, [37j)- Consider the Schrodinger equation = 
E n ty , where V = V- = Vq is a shape invariant potential. If we fix the first level 
of energy Eq = 0, then the excited spectrum and the wave functions are given 
respectively by 

n+l n 

£„ = ^i?K), ¥-\x;a 1 ) = Y[A+(x;a k )^ ) (x;a n+1 ). (1.39) 

k=2 k=l 

The complete statement and proof of theorem ll.2.13l can be found in [37ll31j. 
As consequences we have the following facts (see [37l I3T]): 

tfWfcao) = *<-\x;ax), *^\(x; a ) = A*(x; ao)*^ _) (s; ai). 

71 

=-d 2 x + V-(x;a n ) + Y / RM 

k=l 

has the same spectrum of if for n > 0, where 

H<®=H-, HW=H+, a n = f n (a ). 



n 



H 



(n+l) 



-dl + V+(x;a n )+J2 R ( a k) 

k=l 
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Following [24i 131] we present the list of shape invariant potentials given in expres- 
sion (TOO)) . 



Potential V Name 



\co 2 x 2 — tt Harmonic Oscillator 



\uj 2 r 2 + - (I + |) uj 3D Harmonic Oscillator 



t+i) 



Coulomb 



A 2 + B 2 e- 2ax -2B(A+ §) e- a:E Morse 

-4 2 + 4t - 25 coth ar + A A .^ a £cfcari 

smn ar 

A 2 + + 25 tanh ax - A A X a Rosen-Morse Hyp. 

A 1 cosh' 1 ax 

-A 2 + -C + 2B cot ax + AJ4 2 - Rosen-Morse Trig. 

A 2 + B 2 -A*-Aa + B(2A+a)sinhax g ^ R j 
cosn ax cosh ax J aL 

A 2 + B-+Al+Aa _ B(2A+ a )coBhar g j ff jj 
sinn ar smn ar J Ul 

_ A 2 + B'+A^-Aa _ B(2A-a)sina X g ^ m j 
cos ax cos z ax J a 

_ A 2 + B*+A*-Aa _ B(2A-a) COS ax g ^ m R 
sin ax sin ax J 3 

~(A + B) 2 + ^f^- + Pdschl-Teller 1 

V / cos^ ax sm ax 

(A-B) 2 - + Sgpi. Pdschl- Teller 2 

v 1 cosn ar sinn ar 



(1.40) 



We recall that a good short survey about Darboux transformations can be 
found in [861 . 
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Chapter 2 



Differential Galois Theory 
Approach to Supersymmetric 
Quantum Mechanics 

In this chapter we present our original results of this thesis, which corresponds to 
the Galoisian approach to Supersymmetric quantum mechanics. We start rewrit- 
ing in a Galoisian context some points of the section 11721 chapter [TJ The results 
presented here are also true for any differential field K with field of constants C in 
agreement with definition ll.1.41 We emphasize in the differential fields Kq = C(x) 
and K\ = C(z(x), d x z(x)), where in both cases C = C. 

2.1 Preliminaries 

The main object of our Galoisian analysis is the Schrodinger equation (|1.22[) . which 
now is written as 

£ A := = A*, H = -dl + V(x), V e K, (2.1) 

where K is a differential field (with C as field of constants). We are interested in 
the integrability of equation (|2.ip in agreement with definition 11.1.81 

We introduce the following notations. 

• Denote by A C C the set of eigenvalues A such that equation (|2.ip is integrable 
according with definition 11.1.81 

• Denote by A + the set {A G A n R : A > 0} and by A_ the set {A G A n R : 
A < 0}. 
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• Denote by L\ the Picard-Vessiot extension of C\. Thus, the differential Galois 
group of C\ is denoted by DGal(L\/K). 

The set A will be called the algebraic spectrum (or alternatively the Liouvil- 
lian spectral set) of H. We remark that A can be 0, i.e., DG&\(L\/ K) = SL(2,C) 
VA G C. On the other hand, by theorem[TXgj if A G A then (DGal(L Xo /K)) a C B. 

Definition 2.1.1 (Algebraically Solvable and Quasi-Solvable Potentials). We say 
that the potential V(x) € K is: 

• an algebraically solvable potential when A is an infinite set, or 

• an algebraically quasi- solvable potential when A is a non-empty finite set, or 

• an algebraically non-solvable potential when A = 0. 

When Card(A) = 1, we say that V(x) G K is a trivial algebraically quasi-solvable 
potential. 

Examples. Assume K = C(x). 

1. If V(x) = x, then A = 0, V(x) is algebraically non-solvable, see [551 155] , 

2. If V(x) = 0, then A = C, i.e., V{x) is algebraically solvable. Furthermore, 

DGal(i /A') = e, VG&\{L X /K) =G m , A ^ 0. 

3. If V(x) = x + 5' ^en A = {n : 77, e Z}, V(x) is algebraically solvable. This 
example corresponds to the Weber's equation, see subsection If .1.41 

4. If V{x) = x A — 2.t, then A = {0}, V{x) is algebraically quasi-solvable (trivial). 

Remark 2.1.2. We can obtain algebraically solvable and quasi- solvable potentials 
in the following ways. 

• Giving the potential V, we try to solve the differential equation d^o = V^>o 
expecting to obtain the superpotential W = d x ln^o)- If the superpotential 
exists in the Liouvillian class (the differential equation is integrable), then 
we search A into the Schrbdinger equation = A^&. With this method the 
algebraic spectrum can be the empty set. We will illustrate this in section 
2.2 and 2.3. 

• Giving the superpotential W, we construct the potential V = d x W + W 2 . 
After we can search the algebraic spectrum A in = A^. With this method 
we have that G A, so at least we obtain trivial algebraically quasi-solvable 
potentials. We will illustrate this in section 2.2 and 2.3. 
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• Using integrable parameterized differential equations which can be trans- 
formed into Schrodingcr equations. In this case the parameter of the differen- 
tial equation should coincide with the eigenvalues of H. With this method we 
know previously the algebraic spectrum A, thus, we can obtain algebraically 
solvable or algebraically quasi-solvable potentials depending on A. We will 
illustrate this in section 2.3. 

We are interested in the spectrum (analytic spectrum) of the algebraically solv- 
able and quasi-solvable potentials, that is. spcc(iJ) n A / 0. For example, the 
potential V(x) = \x\ has point spectrum (see [B]) although V(x) is algebraically 
non-solvable. Thus, when spec(-ff ) fl A is an infinite set, in the usual physical ter- 
minology these potentials are called solvable ( or exactly solvable ) potentials, see 
Natanzon 69J. In analogous way, when spec(if) n A is a finite set, the usual def- 
inition in physics of these potentials is quasi- exactly solvable (or quasi-solvable) 
potentials (Turbiner |102j . Bender & Dunne pj], Bender & Boettcher [10], Saad 
et al. [57], Gibbons & Vesselov [3SJ). 

Definition 2.1.3. Let be C, C, pairs of linear differential equations defined over 
differential fields K and K respectively, with^Picard-Vessiot extensions L and L. 
Let if be the transformation such that £ i— ► £, K i— > K and Li->I, we say that: 

1. ip is an iso-Galoisian transformation if 

DGal(A/A') = DGal(Z/AT). 

If L = L and K = K, we say that ip is a strong iso-Galoisian transformation. 

2. ip is a virtually iso-Galoisian transformation if 

(DGal(L/A'))° = (DGal(Z/A'))°. 

Remark 2.1 A. The Eigenrings of two operators £ and £ are preserved under iso- 
Galoisian transformations. 

Proposition 2.1.5. Consider the differential equations 

C := dly + ad x y + by = 0, £:=5gC = rC, a, b, r e K. 

Let k e Q, / e K, a = 2nd x (hi f) and ip be the transformation such that C i— ► C 
The following statements holds: 

1. p is a strong isogaloisian transformation for (t£Z. 

2. p is a virtually strong isogaloisian transformation for k £ Q \ Z. 
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Proof. Assume that B = {y 1,2/2} is a basis of solutions and L is the Picard-Vessiot 
extension of £, B' = {Ci,C2} is a basis of solutions and L is the Picard-Vessiot 
extension of £. With the change of dependent variable y = (e~^ a we obtain 
r = a 2 /4 + d x a/2 — b and for instance K = K. Thus, the relationship between L 
and L depends on a: 

1. If k = n € Z, then 6' = {.f n yi, f n y2} which means that L = L and ip is 
strong isogaloisian. 

2. If /c = ^, with gcd(n, m) = 1, — ^ Z, then 23' = Z™^} which means 
that L is either an algebraic extension of degree at most m of L, and ip is 
virtually strong isogaloisian. or L = L when fm £ ii" which means that ip is 
strong isogaloisian. 

□ 

Remark 2.1.6. The transformation tp in proposition 12.1.51 is not injective, there 
are a lot of differential equations £ that are transformed in the same differential 
equation C. 

As immediate consequence of the previous proposition we have the following 
corollary. 

Corollary 2.1.7 (Sturm-Liouville). Let C be the differential equation 

d x {ad x y) = {Xb - /j,)y, a, be K 7 A,//gC 

in where L, L, C and ip are given as in proposition \2.1.5\ Then either L is a 
quadratic extension of L which means that ip is virtually strong isogaloisian or 
L = L when a? G K which means that ip is strong isogaloisian. 

2.2 Supersymmetric Quantum Mechanics with Ratio- 
nal Potentials 

Along this section we consider as differential field K = <C(x). 
2.2.1 Polynomial Potentials 

We start considering the Schrodingcr equation (|2.1[) with polynomial potentials, 
i.e., V S C[x], see [THl HO]- For simplicity and without lost of generality, we con- 
sider monic polynomials due to the reduced second order linear differential equa- 
tion with polynomial coefficient c n x n + . . . + C\x + Cq can be transformed into 
the reduced second order linear differential equation with polynomial coefficient 
x n + . . . + q\x + go through the change of variable x 1— > n+ ^/c^x. 
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When V is a polynomial of odd degree, is well known that the differential 
Galois group of the Schrodingcr equation (|2.ip is SL(2,C), see [5"o] . 

We present here the complete result for the Schrodingcr equation (|2.1j) with 
non-constant polynomial potential fTheorem l2.2.2p . see also [SJ §2]. The following 
lemma is useful for our purposes. 

Lemma 2.2.1 (Completing Squares, [5]). Every even degree monic polynomial of 
can be written in one only way completing squares, that is, 

2n-l I n-l \ 2 n-l 

Q2n(x) = x 2n + VkX k = x n + 5> fe :r fc +Y, b kx\ (2-2) 

fc=0 V k=0 / fc=0 

where 

_ qin-\ _ 92/1-2 ~ Q n-1 _ g2n-3 ~ 2a K -iQ„_2 

^n — 1 ~ j ei> n —2 ~ , CL n — 3 — , • • • , 



<Zn — 2aia„_i — 2a2a n _2 — ■ ■ ■ , 2 i n 

a = , b Q = q -a , &i = qi - 2a Oi, 



t>n-i — ?n-i — 2aoa„_i — 2aia„_ 



2 



Proo/. See lemma 2.4 in [5J p. 275]. 

□ 

We remark that V(x) as in equation (|2.2|) can be written in terms of the 
supcrpotential W(x), i.e., V(x) = W 2 (x) — d x W(x), when 

n— 1 n— 1 

nx n ~ l + Y] kdkx k ~ l = - \" 6 fc x fc 

fe=l fc=0 

and W(x) is given by 

n-l 

x n + V" a fc x fc . 

fc=0 

The following theorem also can be found in [3 §2], see also [4]. Here we 
present a quantum mechanics adapted version. 



Theorem 2.2.2 (Polynomial potentials and Galois groups, [5]). Let us consider the 
Schrodinger equation (|2.ip . with V(x) <E C[x] a polynomial of degree k > 0. Then, 
its differential Galois group DGal(L\/ K) falls in one of the following cases: 
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1. DGal(L x /K) = SL(2,C), 

2. BGal(L x /K) = E, 

and the Eigenring of H — A is trivial, i.e., £{H — A) = Vect(l). Furthermore, 
DGsl(Lx/ K) = B if and only if the following conditions hold: 

1. V(x) — A is a polynomial of degree k = In writing in the form of equation 

(El. 



2. b n —i — n or — 6„_i — n is a positive even number 2m, m £ Z+. 

3. There exist a monic polynomial P m of degree m, satisfying 

I n-1 \ / n-2 n-1 \ 

dlP m +2 I x n + ^2 a kX k d x P m + I na: 71-1 + ^(fc + l)a k+1 x k - ^ b k x k P m = 0, 

V fc=0 / V fc=0 fc=0 / 



or 



/ n-1 \ / n-2 n-1 \ 

d 2 x P m -2 ix" + a k x k d x P m - nx"- 1 + + l)a k+1 x k + ^ b k x k \ P m = 0. 

V fc=0 / V fc=0 fc=0 / 



In such cases, the only possibilities for eigenf unctions with rational superpotentials 
are given by 

n+1 n-1 fe+1 

*A=P ro e'<*>, or * A = P m e-^), wfere /(*) = + V ^f— . 

n+1 f—f fe + 1 



fe=o 



Proof. See theorem 2.5 in [5j p. 276]. 

□ 



An easy consequence of the above theorem is the following. 

Corollary 2.2.3. Assume that V(x) is an algebraically solvable polynomial poten- 
tial. Then V(x) is of degree 2. 



Proof. Writing V(x) — A in the form of equation (|2.2p we see that 6 n _i — n = 2m 
or — b n -i — n = 2m, where m € Z+. Thus, the integrability of the Schrodinger 
equation with Card(A) > 1 is obtained when b n -i is constant, so n = 1. □ 
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Remark 2.2 A. Given a polynomial potential V{x) such that spec p (H)D A ^ 0, we 
can obtain bound states and normalized wave functions if and only if the potential 
V(x) is a polynomial of degree An + 2. Furthermore, one integrability condition 
of = for A G A is that b 2n must be an odd integer. In particular, if the 
potential 

V(x) = x An+2n + fix 2n , n>0 

is a quasi-exactly solvable, then /i is an odd integer. For this kind of potentials, 
we obtain bound states only when /i is a negative odd integer. 

On another hand, the non-constant polynomial potentials V(x) of degree 4n 
are associated to non-hcrmitian Hamiltonians and VT invariance which are not 
considered here, see [10]. Furthermore, one integrability condition of i/\E f = 
for A G A is that &2n-i must be an even integer. In particular, if the Schrodingcr 
equation 

F* = A*, V(x) =x 4n + nx 2n -\ AG A 
is integrable, then \i is an even integer. 

We present the following examples to illustrate the previous theorem and 
remark. 

Weber's Equation and Harmonic Oscillator. The Schrodingcr equation with poten- 
tial V(x) = x 2 + qix + qo corresponds to the Rehm's form of the Weber's equation 
(|1.19p . which has been studied in section H .1.41 By lemma |2". 2. II we have 

V(x) - A = (x + a ) 2 + b , 00 = 91/2, 6 = qo — 3i/4— A. 

So that we obtain ±60 — 1 = 2m, where m G Z+. If 60 is an odd integer, then 

DGal(i A /A') = B, S(H-X) = Vcct(l), A G A = {±(2m + l)+q -<?i/4 : m G Z+} 

and the set of eigenfunctions is either 

*A = P m e^ x ' 2+qix \ or, * A = P m e~^ x2+qix) . 

In the second case we have bound state wave function and spcc p (i7) fl A = 
spec p (-ff) = {E m = 2m + 1 + q — qf/4 : m G Z + }, which is infinite. The polyno- 
mials P m are related with the Hermite polynomials H m , [2"51 [STIITO] . 

In particular we have the harmonic oscillator potential, which is given in the 
list (|1.40p and where H^> = E^>. Through the change of independent variable 

x 1 — > \f^ x we obtain V{x) = x 2 — 1 and A = ^E, that is, <?i = and go = — !• I n 
this way A = {±(2m + 1) — 1 : m G Z + } and the set of eigenfunctions is either 

* A = P m e* x \ or, * A = P m e-i x \ 
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where as below, DGal(L\/ K) = B and 8 (H — A) = {1} for all A £ A. In the second 
case we have bound state wave function, spec p (H) n A = spec p (-ff) = A + = {2m : 
m £ Z + } and P m = H m . The wave functions of ff\E f = for the harmonic 
oscillator potential are given by 




= mu>. 



Quartic and Sextic Anharmonic Oscillator. The Schrodinger equation with poten- 
tial V(x) = x 4 + q3X 3 + qix 2 + q\x + qo can be obtained through transformations of 
confluent Hcun's equation, which is not considered here. By lemma [2.2.11 we have 

V(x) — A = (x 2 + ol\x + a ) 2 + b\x + bo, 

where ai = (73/2, a = (72/2 — &i/2, &i = q\ — 2a^ai and b a = q — Oq — A. So that 
we obtain ±61 — 2 = 2m, where m g Z + . If A 7^ 0, then b\ is an even integer, 
P m satisfy the relation ((OJ) and T)G&\{L\/K) = B for all HA. The set of 
eigenfunctions is cither 

* A = P m e-+— +aoX , or, * A = P m e \ 3 J, 

where A and m are related, which means that A is finite, i.e., the potential is 
algebraically quasi-solvable. In particular for (73 = 2il, q^ = I 2 — 2k, qi = 2t(lk — J) 
and qo = 0, we have the quartic anharmonic oscillator potential, which can be 
found in [TU]. 

Now, considering the potentials V(x,\i) — x 4 + 4x 3 + 2x 2 — fix, again by 
lemma [2.2. II we have that 

V(x, fi) - A = {x 2 + 2x - l) 2 + (4 - fi)x - 1 - A, 

so that ±(4 — fx) — 2 = 2n, where neZ + and in consequence (1 £ 2Z. Such fi can 
be either fi = 2 — 2n or fi = 2n + 6, where n £ Z+. By theorem 12.2.21 there exists 
a monic polynomial P n satisfying respectively 

d 2 P n + {2x 2 +4x-2)d x P n + ((n-2)x + 3 + X)P n = 0, /x = 2 - 2n, or 

d 2 x P n - (2x 2 +Ax- 2)d x P n + ((/J - 6)x - 1 + A)P„ = 0, /j = 2n + 6 

for A ^ 0. This algebraic relation between the coefficients of P n , fi and A give us 
the set A in the following way: 

1. Write P n = x n + c„_ix Il_1 + . . . + cq, where Ci are unknown. 

2. Pick fi and replace P n in the algebraic relation (|1.6[) to obtain a polynomial 
of degree n with n+ 1 undetermined coefficients involving Co, . . . c rl _i and A. 
Each of such coefficients must be zero. 
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The term n + 1 is linear in A and c„_i, thus we write c„_i in terms of A. 
After of the elimination of the term n + 1, we replace Cn—x in the term n 
to obtain a quadratic polynomial in A and so on until arrive to the constant 
term which is a polynomial of degree n + 1 in A (Q„+i(A)). In this way, 
A = {A : Qti+i(A) = 0} and Co, ■ ■ ■ , Cn-i are determined for each value of A. 

For fi = 2n + 6, we have: 

n = 0, V(x,6), P = l, A = {1} 

n=l, V(x,S), Pi = x + 1 =F V2, A = {3±2V2} 



and the set of eigenfunctions is 

In the same way, we can obtain A, P n and ^a,^ for fj, = 2 — 2n. However, 
we have not bound states, spec p (-ff) n A = 0, DGal(iA/^') = B and £(H — A) = 
Vect(l) for all A g A. 

The well known sextic anharmonic oscillator x e + q§x° + ■ ■ ■ + qix + qo can be 
treated in a similar way, obtaining bound states wave functions and the Bcndcr- 
Dunnc orthogonal polynomials, which corresponds to Q n+ i(A), i.e., we can have 
the same results of [TTJ [351 ISZ] ■ The Schrodingcr equation with this potential, un- 
der suitable transformations, also falls in a confluent Hcun's equation. 



2.2.2 Rational Potentials and Kovacic's Algorithm 

In this subsection we apply Kovacic's algorithm to solve the Schrodinger equation 
with rational potentials listed in the previous chapter (equation (|1.40[) ). 

Three dimensional harmonic oscillator potential: 

we can see that Schrodinger equation (equation (|1.22p ) for this case can be written 

as 

^((H'+^-H)-*)-- 

By the change I » / — ) r we obtain the Schrodinger equation 



^ = ( r a + £(£+l)_(2^ + 3)-A > )* 1 X = -E 
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and in order to apply Kovacic's algorithm, we denote: 

2 , 



R = r z 



(2£ + 3) - A. 



We can see that this equation could fall in case 1, in case 2 or in case 4 
(of Kovacic's algorithm). We start discarding the case 2 because by step 1 (of 
Kovacic's algorithm) we should have conditions C2 and 003, in this way we should 
have E c = {2,4 + 4£, —A£} and Eoo = {—2}, and by step 2, we should have that 
n = —4 ^ Z + , so that D = 0, that is, this Schrodinger equation never falls in case 
2. Now, we only work with case 1; by step 1, conditions ci and 003 are satisfied, 
so that 



Vr 



0, a 



± _1±(2£+1) 



± t(A + 2£ + 3) - 1 



By step 2 we have the following possibilities for n G Z + and for A £ A: 



A++) 


n = 


«i - 


a+ 


= -i (4^ + 6 + A) 


, A = 


-2n-4£-6, 


A+-) 


n = 


a+ - 




= -|(4 + A), 


A = 


-2n-4, 


A-+) 


n = 




at 


A 
2 ' 


A = 


2n, 


A-) 


n = 




a o 


= i (4£ + 2 + A) , 


A = 


2n - 4g - 2, 


where A ++ U A + 


_UA_ 


-+U A 




= A, which means 


that A 


= 2m, m G Z. Now, 



A G A, the rational function cu in Kovacic's algorithm is given by: 



A++) 


u = r + 


t+i 




= r 2 H 




+ {2£ + 


3) 


+ 2n 


r 1 






A+-) 


w = r — 


£ 
r * 


-Rn 


= r 2 - 


1(1+1) 


~ {2£~ 


1) 


+ 2n 


A-+) 


w = — r 






= r 2 - 


1(1+1) 


- (2£ + 


3) 


- 2n 








A—) 


uj = —r 


r ' 


-Rn 


= r 2 ^ 


£(£+1) 

r 2 


+ {2£- 


1) 


- 2n 



where R n is the coefficient of the differential equation C n := d 2 ^ = R n ^, which 
is integrablc for every n and for every A G A we can see that DGal(L„/ K) = 
DGal(L x /K), where C x := = A* and A G A. 

By step 3, there exists a polynomial of degree n satisfying the relation () 1 .G|) : 

A++) d 2 r P n + 2 (r + £±1) d r P n - 2nP n = 0, A G A_, 

A+_) d 2 r P n + 2 (r - i) 3 r P n - 2nP n = 0, A G A_, 

A_+) d 2 P n + 2(-r + ^)d r P n + 2nP n = 0, A G A+, 

A__) d 2 P n + 2(-r-^)d r P n + 2nP n = 0, A G A. 



2.2. Supcrsymmctric Quantum Mechanics with Rational Potentials 



53 



These polynomials exists for all A G A when their degrees are n 6 2Z, while 
for n G 2Z + 1, they exists only for the cases A |_) and A ) with special con- 
ditions. In this way, we have obtained the algebraic spectrum A = 2Z, where 
A ++ = 4Z_, A + _ = 2Z_, A_ + = 4Z+, A = 2Z. 

The possibilities for eigenfunctions, considering only A G 4Z, are given by 



A++) 


*n(r) 


= r< 


' +1 P 2 n(r)e^, 


A G A_ 


A+-) 


*n(r) 


= r 




A G A_ 


A-+) 


*n(r) 


= r' 


' +1 P 2)l (r)e^, 


A G A+ 


A-) 


*n(r) 


= r 




A G A. 



To obtain the point spectrum, we look ^„ satisfying the bound state conditions 

(equation (|1.23[) ) which is in only true for A G A |_. With the change r i— > •v/fr, 

the point spectrum and ground state of the Schrodinger equation with the 3D- 
harmonic oscillator potential are respectively spec p (i?) = {E n : n G Z + }, where 
E n = 2nu, being w the angular velocity, and 

— \ <+i 



d-r 

The bound state wave functions are obtained as ^ n = P 2 n^o- Now, we 
can see that DGa\(L /K) = B and £{H) = Vect(l). Since = P 2 „*o, for all 
A G A we have that DGal{L\/K) = B and £{H - A) = Vect(l). In particular, 
T>G&\(L\/K) = B and £(H - A) = Vect(l) for all A G spec p (£T), where A = %E. 

We remark that the Schrodinger equation with the 3D-harmonic oscillator 
potential, through the changes r i— > |wr 2 and ^ i— > v 7 ^) m a Whittakcr 
differential equation (equation (|1.16p ) in where the parameters are given by 

(2£ + 3)u + 2E 1 1 

k = - , u=-£+-. 

Applying theorem ll.l.31i we can see that for integrability, ±k± /i must be a half 
integer. These conditions coincides with our four sets A±±. 

Coulomb potential: 

nr) = _!! + *±i) + _/_, <ez . 

w r r 2 4(£+ 1) 2 

we can see that the Schrodinger equation (equation (|1.22[) ) for this case can be 
written as 



(t(t+X\ e 



2 e 4 



54 



Chapter 2. Differential Galois Theory Approach to SUSY QM 



By the change r \— > 2 ^t"^ r we obtain the Schrodinger equation 

( 

and in order to apply Kovacic algorithm, we denote 



Firstly we analyze the case for A = 1: we can see that this equation only 
could fall in case 2 or in case 4 of Kovacic's algorithm. We start discarding the 
case 2 because by step 1 we should have conditions C2 and 003. In this way we 
should have E c — {2, 4 + 41, — 4€} and = {1} and by step 2, we should have 
that n ^ Z. Thus, n Z+ and D = 0, that is, the differential Galois group of this 
Schrodinger equation for A — 1 is SL(2,C). 

Now, we analyze the case for A 7^ 1: we can see that this equation could fall 
in case 1, in case 2 or in case 4. We start discarding the case 2 because by step 1 
we should have conditions C2 and 003, so that we should have E c = {2,4 + 4^, —M} 
and Eqq = {0}. By step 2, we should have that n = 2£ £ Z+, so that D = {2£} 
and the rational function 8 is 8 = =^ , but we discard this case because only could 
exists one polynomial of degree 11 for a fixed £, and for instance, only could exist 
one eigenstate and one eigenfunction for the Schrodinger equation. 

Now, we only work with case 1, by step 1, conditions C2 and 003 are satisfied. 
Thus, 



Q ± = I±1^±±) 



VI -A, a*=T- 



By step 2 we have the following possibilities for n £ Z + and for A 6 A: 



A++) n = a+ -a£ = -(l + 1) (l + 7 ^=) 



A = 1 



, 2 



l+l+n 
2 



2 

-1 



A -+) n = a x -a+ = (e+l)(-^=-l), A=l L+ „ 
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We can see that A € A_ when A < 0, while A e A + when < A < 1. Furthermore: 



A_, I < 



A++) t<-l, Ae 



A+_) £>0, Ae 



A_+) £eZ, Ae 



A__) e>0, Ae 



-n — 2 
2 



a+, =^<e<-i 



n-l 



A_, l> 2 



A+, 0<€<^± 



A_, I < 



-n-2 



2 

A+, l>-\ 



A_, 



A+, 0<^<^ 



In this way, the possible algebraic spectrum can be A = A++ U A^ U A |_ U A_ 

that is 



K =^-\iTTT-n :nez+ r v-\j^ :neZ+} > (2 - 3) 



Now, for A e A, the rational function uj is given by: 



a ++ ) -=^+^, ^a ++ , R^^-^ + ^y, 

A G Ah — , i? n = £(^2(£±I> + (f±l) 2 , 



A_ + ) w = _*£_ + *+I A e A_+, R n = e -^± -Hi£±ii + 

2 



- r 2 


r 


e(t+i) 

- r'± 


1.(1+1) 

r 


1(1+1) 


i(i+i) 

r 


1(1 + 1) 


i(t+i) 

r 



where R n is the coefheient of the differential equation = R n ^, which is in- 
tegrable for every n. 
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By step 3, there exists a polynomial of degree n satisfying the relation (|1.6p . 
A ++ ) ^» + 2(^ + ^)^+2^(1 + ^)^ - 0. 



^ + 2(|M_l) 9r P„ + ?(ai)(l-|±i 

A- + ) d % r Pn + 2 {-M~ n + ^) d r P n + &±H (l 



0, 



i+l+n J " 

A-) ^P„ + 2 (-|±i-i)9 r P„ + ^±i)(l + |±i)p„ = 0. 

These polynomials exists for every A £ A when n G Z, but Pq = 1 is satisfied 

only for A £ A |_. In this way, we have confirmed that the algebraic spectrum A 

is given by equation (|2.3p . 



The possibilities for cigenfunctions arc given by 

A + + ) *„(r) = r*+ 1 P n (r)/„(r)e r , / n (r) = e«+!+« , Aej 

n+1 f 

A + _) *„(r) = r- e P n (r)f n (r)e r , fn(r) = , \ £ i 

A_ + ) * n (r) = r^+ 1 rv,(r)/„(r)e- r , /„ (r) = e , Aej 



n+1, 



{ 



A- , 


I 


< 


-n— 2 
2 


A + , 




2 


~ 2 < e < 


A_ , 


I 


> 


2 


A + , 





< 


* - ~T~ 


A_ , 


I 


< 


-n—2 
2 


A + , 


I 


> 


- 1 


A_ , 


I 


> 


2 


A + , 





< 


« < 2=1 
- 2 



but should satisfy the bound state conditions (equation (|1.23[) ) which is only 
true for A e A h n A + , so that we choose A |_ n A + = spec p (H), that is 

SpeC ^)-{l-(^r) 2 : nSZ,, <>-l}. 

2 

By the change r 1— > 2 (i+i) r ' ^ ne point spectrum and ground state of the Schrodingcr 
equation with Coulomb potential are respectively 



spec p (if) = {E n :neZ + }, E„ 



1 



4 \{£+l) 2 (e+l + n) 2 
and 



e+i 
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The eigenstates are given by ^ n = P n fn^o, where 

Tie 2 r 

f n ( r ) = e 2(f + l + r,)(f+l) _ 

Now, we can see that DGal(L /K) = B and£(H) = Vcct(l). Since = F 2 n/n*o, 
for all A e A we have that BGal{L x /K) = 1 and £(H - A) = Vect(l). In par- 
ticular, DGal(i A /A') = B and £(H - A) = Vcct(l) for all £ € spec p (£T), where 

^ = 4(f+l) 2 ^- 

We remark that, as in the three dimensional harmonic oscillator, the Schrodingcr 
equation with the Coulomb potential, through the change 

a/-4 (£+l) 2 £ + e 4 

r i— > — r, 

^ + 1 

falls in a Whittakcr differential equation (equation Ql.lfip ) in where the parameters 
are given by 

-4(£+l) 2 £ + e 4 ^ 

Applying theorem 1 1.1.31) we can impose ±k±/i half integer, to coincides with our 
four sets A±±. 

Remark 2.2.5. By direct application of Kovacic's algorithm we have: 

• The Schrodinger equation (|2.1|) with potential 

V(x) = ax 2 + 4r 
ar 

is intcgrablc for A € A when 

-a = 0, 6 = /x(/i+l), ,ueC, A = C, 
- a = 1, & = 0, A e A = 2Z+ 1, 
-a=l, 6 = ^ + 1), feZ*,A = 2Z + l. 

• The only rational potentials (up to transformations) in which the elements 
of the algebraic spectrum are placed at the same distance, belongs to the 
family of potentials given by 



k— — oo 



In particular, the set A for the harmonic oscillator (a = 1, b = 0) and 3D 
harmonic oscillator (a — 1, b — £(£ +1)) satisfies this. 



58 



Chapter 2. Differential Galois Theory Approach to SUSY QM 



Proposition 2.2.6. Let C\ be the Schrddinger equation (|2.1|) with K = C(x) and 
Picard-Vessiot extension Lx- IfT)Gal(L /K) is finite primitive, thenDG&\{L\/ K) 
is not finite primitive for all A £ A \ {0}. 

Proof. Pick V £ C(x) such that £q falls in case 3 of Kovacic algorithm, then 
olio, > 2. Assume t, s £ C[x] such that V = f, then deg(i) > deg(s) + 2 and 
V — A = a ~ Af . Now, for A ^ we have that deg(s — At) = deg(t) and therefore 
o(V — A)oo = 0. So that for A ^ 0, the equation C\ does not falls in case 3 of 
Kovacic algorithm and therefore DGal(L\/K) is not finite primitive. □ 

Corollary 2.2.7. Let C\ be the Schrddinger equation (|2.1|) with K = C(x) and 
Picard-Vessiot extension L\. Lf Card(A) > 1, then there is either zero or one 
value of X such that T)Gal{L\/ K) is a finite primitive group. 

Proof. Assume that Card(A) > 1. Thus, by proposition 12.2.61 the Schrodingcr 
equation does not falls in case 3 of Kovacic's algorithm. □ 



It seems that the study of the differential Galois groups of the Schrodingcr 
equation with the Coulomb potential has been analyzed by Jean-Pierre Ramis 
using his summability theory in the eighties of the past century, see [82] . 



2.2.3 Darboux Transformations 

Here we present a Galoisian approach to Darboux transformation, Crum iteration 
and shape invariant potentials. We denote by W(yi, . . . ,y n ) the Wronskian 



i,---,y n ) 



0':, 



yx 



tin 



u x yn 



by DT the Darboux transformation, by DT„ the n iteration of DT and by CI„ 
the Crum iteration. Also we use the notation of subsection 11.21 We recall that 
K = C(x) and in case of other differential fields, as usually is considered along 
this thesis, we mean the smallest differential containing the coefficients of the 
linear differential equations. 

Theorem 2.2.8 (Galoisian version of DT). Assume H± — &l + V±(x) and A^O. 
Let £\ given by the Schrddinger equation H_^?(' = A^' - ^ with V-(x) £ K and 
C\ given by the Schrddinger equation H + ^^ +s> = A4' ( - + - ) with V+(x) £ K. Let 



DT be the transformation such that C 
following statements holds: 

■[<.-) a2 



C, V- h-> v+, * ( 



#(+). Then the 



i) DT(^_) = V + = *\~>d 2 x I ) + Ax = V- - 2d 2 x (\n*{- ) ) 



DT(5 , ^ i ') = = ^pj, where ty^J is a particular solution ofC\ t , \\ £ A. 



,(+) 



-,(-) 
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ii) DT(*<->) = *{ +) = a a *<-> - fl a (ln *<->)*<-> = \ A ^ Ai, to/iere 

' is the general solution of C\ for A £ A \ {Ai} and is the general 

solution of C\ also for A G A \ {Ai}. 

In agreement with the previous theorem we obtain the following results. 

Proposition 2.2.9. DT is isogaloisian and virtually strong isogaloisian. Further- 
more, if d x (ln^ Xi ) G K, then DT is strong isogaloisian. 

Proof. Let K, L\ be the differential field and the Picard-Vessiot extension of the 
equation C\. Let K, L\ the differential field and the Picard-Vessiot extension 
of the equation £ A . Due to DT(F_) = V + = 2W 2 - V- - 2A X , where W = 
-^(ln*^'), we have K = K (d x {m ^J)) ■ B ^ theorem [TXTU1 we have that 
the Riccati equation d x W = V- — W 2 has one algebraic solution, in this case 
W = — 9a; (In Let (^L be a basis of solutions for equation C\ and 

(^(ia)> ^(^i)) a basis °f solutions for equation C\. Since the differential field for 
equation £ A is K = K(d x (]n i b^~ i ')}, we have that L = K {^c[\y and 

Z = *w } > = )f *w a a (in*i7>)> 

for A = Ai and for A ^ \\. Since dxiln^^) l& algebraic over K, then 

(DGal(L A /A'))° = (DGal(Z A /A'))°, DGal(L A /AT) = DGal(Z A /A'), 

which means that DT is a virtually strong and isogalosian transformation. 

In the case d x (\n^ x ^) £ K, then K = K and L = L, which means that DT is a 

strong isogalosian transformation. □ 

Proposition 2.2.10. Consider £ A := f/_ — A and £ A := H + — A such that DT(i/_ — 
A) = H + — A. The eigenrings of £ A and £ A are isomorphic. 

Proof. Assume £(£ A ) and £(£ A ) the eigenrings of £ A and £ A respectively. By 
proposition 12.2.91 the connected identity component of the Galois group is pre- 
served by Darboux transformation and for instance the eigenrings is preserved by 
Darboux transformation. Now, suppose that T £ £(£ A ), Sol(£ A ) and Sol(£ A ) the 
solutions space for = and £ A ^(~) = respectively. To transform £ (£ A ) 

into £ (£ A ) we follows the diagram: 
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Sol(£ A ) 
Sol(£ A ) 



Sol(£ A ) 
Sol(£ A ) 



T = ATA^raod(£\) G £ (£ A ), 



where A^ and A are the raising and lowering operators. 



□ 



Example. Assume the Schrodinger equation C\ with potential V = V- = 0, which 
means that A = C. If we choose Ai = and as particular solution ' = x, then 
for A ^ the general solution is given by 



c 2 e 



Applying the Darboux transformation DT, we have DT(£ A ) = C\, where 

DT(^_) = V+ = % 



and for A ^ 



^(+) = c^V^Xx- l)e^ 
A x 



- A.r 



-\x 



We can see that K = K = C(x) for all A G A and the Picard-Vessiot extensions 



are given by Lq — Lq = C(x), L\ = L\ 



for A G C* . In this way, we 



have that DGal(L /A) = BGal{L /K) = e; for A ^ 0, we have DGal(L A /A') = 
DGal(L A /A') = G m . The cigenrings of the operators £o and £ A are given by 

£(£o) = Vcct (l, xd x , xd x — 1, x 2 d x — x) , 



l,x*d x -2x 3 ,^ + ^ 



1 

i 



and for A 7^ 



£(£ A ) = Vect (1, d x ) , £(£ A ) = Vect I 1, - I A + ) d x - \ 



1 



where £ A := £ A *(") = and £ A := Z x ^ {+) = 0. 

Applying iteratively the Darboux transformation, theorem 12.2.81 and by 
propositions 12.2.91 12. 2. 101 we have the following results. 
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Proposition 2.2.11 (Galoisian version of DT„). Let be A ^ 0, £\ n ' given by 
#(n)^(n) = A*("), V n € JT„, = K, V = V-, ff (0) = H_, *(°) = and 
the Picard-Vessiot extension of . Let £^ n+1 ^ given by if( n + 1 )\J/( n + 1 ) = 
A*^ 1 ), y„ +1 e A" n+ i. Lei DT„ suc/i £/ia£ £ { ™ ] i-> £^ +1 \ F„ >-> K+i, ^ ^ 
and L^" +1 ' ) i/ie Picard-Vessiot extension of C^ +1 \ Then the following 
statements holds: 

i) DT n (y_) = DT(K) = K+i = V n -2dl (ln*W) = V_-2£Lo^ (l**?)' 

where ^ \^ is a particular solution for X = X^, k = 0, . . . , n. In particular, if 
A n = Ao and A = C, then there exists ^ swc/i that V n ^ V n -2, with n > 2. 

ii) DT(*W) = DT„(*l->) = ^ = O^-^-^f = where 

V I / ^ is a general solution for X £ A \ {A n } of C^™K 

iii) if n+1 = K n (d x (\nV<£)y 

iv) DT„ is isogaloisian and virtually strongly isogaloisian. Furthermore, ifd x (ln^^ ) G 

-fsT n i/iera DT„ zs strongly isogaloisian. 

v) TTie eigenrings of — X and — A are isomorphic. 

Proof. By induction on theorcm l2.2.8l we obtain i) and ii). By induction on propo- 
sition 12.2.91 we obtain iii) and iv) . By induction on proposition 12.2.101 we obtain 
v). 

□ 



Example. Starting with V = 0, the following potentials can be obtained using 
Darboux iteration DT„ (see pHIT7] ). 

n(n-l)b 2 . m 2 n(n- l)(b 2 - a 2 ) 



(bx + c) 2 ' 71 (acosh(ma;) + 6sinh(mx)) 2 ' 

jjn y _ -4a6m 2 n(n - 1) IV \y- m 2 n(n ^ l)(b 2 + a 2 ) 



(ae ml + oe~ ma: ) 2 ' (acos(mx) + 6sin(ma;)) 2 

In particular for the rational potential given in /), we have K = K n = <C(x) and 
for A„ = A = 0, we have 

*n n) = „ Cl x + C2(bx + c) n+1 , so that DGal(Lo/^0 = DGal(4 ,l) /A") = e, 
(bx + c) n 



= Vcct (l,xd x - l,x 2n+2 d x - (n+l)x 2n+ \ —„ 



n 



2n 1 x 2n+l 



G2 
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whilst for A 7^ and A„ = 0, the general solution \E , , r ^ is given by 

\I/ A (x) = cif n (x, X)h n (sin(VXx) + c 2 g n (x, A)j„(cos 
where f n ,g n , h n ,j„ G C(x), so that 

DGal(L A /A") = DGal(L^ ri) /A") = G m , 

and 

dim c £(# - A) = dim c £ (# (n) - A) = 2. 
Proposition 2.2.12 (Galoisian version of CI n ). Consider C\ given by = X^>, 
H = -dl + V, V G K, such that Card(A) > n for a fixed n G Z+. Lei £ A " } 6e 
gwen 6?/ = A*W, w/iere ifM = 9| + ^ y n e irf ci„ 6e the 

transformation such that C\ 1— > £ A > V 1— » V^, (f ^, . . . , ^A n i ^a) |—> ^x > where 
for k = 1, . . . , n and the equation C\, the function ^ \ is the general solution for 
A 7^ Afe and ^ \ k is a particular solution for A = A^. Then the following statements 
holds: 

i) CI„(£ A ) = 4"' where Gl n (V) =V n = V - 2fig (lnW(* Al> • ■ ■ , #aJ) and 



CI n (* A ) = *™ = 



_ , Tf (n) _ W^ Xll -..^X n ^x) 



where is i^e general solution of C J . 

ii) K n = K(d x (lnW(y Xl ,...,y Xn ))). 

iii) CI„ is isogaloisian and virtually strongly isogaloisian. Furthermore, if 

d x (]nW(^ x ^---,^x n ))€K n , 

then CI n is strongly isogaloisian. 

iv) The eigenrings of H — A and — A are isomorphic. 

Proof. By induction on theorem 12.2.81 we obtain i) . By induction on proposition 
12.2.91 we obtain ii) and iii). By induction on proposition 12 . 2 . lUl we obtain iv). 

□ 

Example. To illustrate the Crum iteration with rational potentials, we consider 
V = Jf. The general solution of £ A := H^> = A* is 

de kx (kx— 1) C2e~ kx (kx + 1) 2 

1 , A = ~k , 

X X 

the cigenfunctions for Ai = — 1, and A2 = —4, are respectively given by 

e~ x (x + l) e~ 2x (2x + l) 

= -, *_ 4 = -. 

x 2x 
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Thus, we obtain 

Ch(V) = v 2 



(2x + 3) 2 

and the general solution of £ A 2) := H^^^ = A* (2) is 

PT AT/ ^ - *( 2 ) - °1 (k(2x + 3)-2)e k - c 2 (2 + fc(2x + 3))e- fa _ 2 
CI 2 (* A )-* A ^-^ + , A--fc. 

The differential Galois groups and cigenrings are given by: 

DGal(L /X) = DGal(4 2) M") = e, dim c £(.ff) = dim c £(£r< 2 >) = 4, 
and for A 7^ 

DGal(i A /iT) = DGal(L A 2) /^) = <&m, dim c £(-ff - A) = dim c £(H (2) - A) = 2. 

Proposition 2.2.13. The super symmetric partner potentials V± are rational func- 
tions if and only if the superpotential W is a rational function. 

Proof. The supersymmetric partner potentials V± are written as V± = W 2 ±d x W. 
We start considering the superpotential W £ C(x), so trivially we have that V± £ 
C(x). Now assuming that V± £ C(x) we have that d x W £ C(x) and W 2 £ C(x), 
which implies that Wd x W £ C(x) and therefore W £ C(x). 

□ 



Corollary 2.2.14. The superpotential W £ C(x) if and only i/DT is strong isoga- 
loisian. 

Proof. Assume that the superpotential W £ C(x). Thus, by proposition ^. 2. 91 DT 
is strong isogaloisian. Now, assume that DT is strong isogaloisian. Thus, V± £ C(x) 
and by proposition 12 . 2 . 1 31 we have that W £ C(x). 

□ 

The following definition is a partial Galoisian adaptation of the original defi- 
nition given in [37J (K = C(x)). The complete Galoisian adaptation is given when 
K is any differential field. 



Definition 2.2.15 (Rational Shape Invariant Potentials). Assume V±(x; //) £ C(x; //), 
where [i is a family of parameters. The potential V = V- G C(x) is said to be 
rational shape invariant potential with respect to [i and E = E n being n £ Z + , if 
there exists / such that 

n+l 

V + (x;a ) = V_(x;a 1 ) + R(ai), a 1 = f(a ), E n = y^R(a k ), E =0. 

k=2 
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Remark 2.2.16. We propose the following steps to check whether V G C(x) is 
shape invariant. 

Step 1. Introduce parameters in W(x) to obtain W(x; fi), write V±(x; /z) = W 2 (x; /z)± 
d x W(x] /i), and replace [i by ao and a\. 

Step 2. Obtain polynomials V € C[x; ao, oi] and Q G C[x; ao, ai] such that 

a „, / \ T/ / \ \ P(x;a ,ai) 

a x (F + (x;a ) - VL(x;ai)) = — ? -. 

Q(x;a ,ai) 

Step 3. Set 7 ? (x;ao,ai) = 0, as polynomial in x, to obtain ax in function of ao, 
i.e., ai = /(ao). Also obtain R(a\) — V+(x;ao) — VL(x;ai) and verify that 
exists fc G Z + such that i?(ai) + • • ■ + R(a,k) ^ 0. 

Example. Consider the supcrpotcntial of the three dimensional harmonic oscillator 
W(r;£) = r — By step 1, the supersymmetric partner potentials are 

By step 2, we have d r (y+(r;a ) - VL(r;ai)) = -2 a ° +3a °~a'" ai+2 . By step 3, 
(ao + l)(ao + 2) = ai(ai + 1), so that a\ = /(ao) = ao + 1, a n = /(a n _i) = a + n, 
i?(ai) = 2. Thus, we obtain the energy levels E n = 2n and the wave functions 
^ { n\r\i) = A^(r;t)---A^(r;£ + n- l)* _) (r; £ + n), compare with [3T]- 

By theorem [2T81 and propositions [2~2~9l [2XTU1 and I2XT31 we have the fol- 
lowing result. 

Theorem 2.2.17. Consider C„ := H^ 1 *^ = E n ^' with Picard-Vessiot extension 
L n , where n G Z_|_. If V — V— G C(x) is a shape invariant potential with respect 
to E = E n , then 

DGal(L n+1 //v) =DGal(L„/A'), E{H - E n+l ) ~ E{H - E n ), n > 0. 



Remark 2.2.18. The differential automorphisms a commutes with the raising and 
lowering operators A and A^ due to W G C(x). Furthermore the wave functions 
^>n ' can be written as ' = P n fn^o \ where P n is a polynomial of degree n 
in x and f n is a sequence of functions being /o(x) = 1 such as was shown in the 
case of Harmonic oscillators and Coulomb potentials. 
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2.3 The Role of the Algebrization in Supersymmetric 
Quantum Mechanics 

In supersymmetric quantum mechanics, there exists potentials which are not ra- 
tional functions and, for this reason, it is difficult to apply our Galoisian approach 
such as in section [2721 In this section we give a solution to this problem presenting 
some results concerning differential equations with non-rational coefficients. For 
these differential equations it is useful, when is possible, to replace it by a new dif- 
ferential equation over the Riemann sphere P 1 (that is, with rational coefficients). 
To do this, we can use a change of variables. The equation over P 1 is called the 
algebraic form or algebrization of the original equation. 

This algebraic form dates back to the 19th century (Liouville, Darboux), but 
the problem of obtaining the algebraic form (if it exists) of a given differential 
equation is in general not an easy task. Here we develop a new method using the 
concept of Hamiltonian change of variables. This change of variables allow us to 
compute the algebraic form of a large number of differential equations of different 
types. In particular, for second order linear differential equations, we can apply 
Kovacic 's algorithm over the algebraic form to solve the original equation. 

The following definition can be found in [T31 [UJ EE] • 

Definition 2.3.1 (Pullbacks of differential equations). Let £i £ Ki[d z ] and £2 £ 
K2[d x ] be differential operators, the expression £2 65 (d x + v) refers to the operator 
whose solutions are the solutions of £2 multiplied by the solution e~ ■> vdx of d x + v. 

• £2 is a proper pullback of £1 by means of / £ K% if the change of variable 
z = f(x) changes £1 into £2. 

• £2 is a pullback (also known as weak pullback) of £1 by means of / £ K2 
if there exists v £ K 2 such that £2 ® (d x + v) is a proper pullback of £1 by 
means of /. 

In case of compact Riemann surfaces, the geometric mechanism behind the 
algcbrization is a ramified covering of compact Riemann surfaces, see |68[ 167] . 

2.3.1 Second Order Linear Differential Equations 

Some results presented in this subsection also can be found in [5J §2 ]. 

Proposition 2.3.2 (Change of the independent variable, [5]). Let us consider the 
following equation, with coefficients in <C(z): 

C z := dig + a(z)d z y + b(z)y = 0, (2.4) 

and C(z) c — ► L the corresponding Picard-Vessiot extension. Let (K,6) be a differ- 
ential field with C as field of constants. Let 8 £ K be a non-constant element. 
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Then, by the change of variable z = 9(x), equation (|2.4|) is transformed in 
£ x :=d*r+(a(9)d x e-?^Jd x r + b{9)(d x 6) 2 r = 0, 8 X = 5, r = yo9.(2.5) 

Let Kq C K be the smallest differential field containing 9 and C. Then equa- 
tion (|2.5|) is a differential equation with coefficients in Kq. Let Kq Lq be the 
corresponding Picard-Vessiot extension. Assume that 

CO) -> K Q , z^>6 

is an algebraic extension, then 

(DGal(Ao/A ))° = (DGal(L/C(z)))°. 

Proposition 2.3.3. Assume C x and C z as in proposition \2.3.^ Let tp be the trans- 
formation given by 

z i— > 9(x) 

<p : 

d * » ate 6 - 

Then DGal(L /A ) ~ DGal(X/A' n L) C DGal(L/C(z)). Furthermore, ifK nL 
is algebraic over C(z), then (DGal(L /A ))° ~ (DGal(Z/C(>))) . 

Proof. By Proposition 12.3.21 the transformation <p leads us to 

C(z) ~ p(C(«)) A , 

that is, we identify C(z) with y>(C(z)), and so that we can view C(z) as a subficld 
of Ao and then by the Kaplansky's diagram (see [53l 1106] ). 

Lo . 

. DGal(L /if ) 



L Kq 



DGal(L/(Ln^o))^ % 

L^Kq 



C(z) 

so that we have 

DGal(i /A ) ~ DGal(L/A' fll) C DGal(L/C(z)) 
and if Kq n L is algebraic over C(z), then 

(DGal(L /A ))° ~ (DGal(L/C(z)))°. 

□ 
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Along the rest of this section we write z = z(x) instead of 9. 



Remark 2.3 .4 (Hard Algcbrization). The proper pullback from equation (|2.5jl to 
equation (|2.4p is an algcbrization process. Therefore, we can try to algebrize any 
second order linear differential equations with non-rational coefficients (proper 
pullback) if we can put it in the form of equation (|2.5p . To do this, which will be 
called hard algebrization, we use the following steps. 

Step 1. Find (d x z) 2 in the coefficient of y to obtain d x z and z. 

Step 2. Divide by (d x z) 2 to the coefficient of y to obtain b{z) and check whether 
beC(z). 

Step 3. Add (d 2 z)/d x z and divide by d x z to the coefficient of d x y to obtain a{z) 
and check whether a £ C(z). 

To illustrate this method, hard algebrization, we present the following example. 



Example. In [92l p. 256], Singer presents the second order linear differential equa- 
tion 

a;(lna; + 1) 

To algebrize this differential equation we choose (d x z) 2 = (In a; + l) 2 , so that 
d r z = lnx + 1 and for instance 



(In a; + \)dx = xhix, b(z) = — 1. 
Now we find a(z) in the expression 
a(z)(\nx + V 



x(hix + l) xQnx+1) 

obtaining a(z) — 0. So that the new differential equation is given by d 2 y — y = 0, 
in which y{z(x)) = r(x) and one basis of solutions of this differential equation is 
given by (e z ,e~ z ). Thus, the respective basis of solutions of the first differential 
equation is given by (e xlnx ,e~ xlnx ). 

In general, this method is not clear because the quest of z = z(x) in b(z)(d x z) 2 
can be purely a lottery, or simply there is not exists z such that a(z), b(z) £ C(z). 
For example, the equations presented by Singer in [2H p. 257, 261, 270] and given 
by 

„, q=2xln 2 x =F 2xlnx — 1 „ — 2xln 2 x — 3xlna; — x =F 1 

dlr + Z d x r + = 0, 

xmx + x x mx + x 

4.x In x ~t~ 2x 

d 2 r+ d x r— r = 0, (x 2 In 2 x)d 2 r+{xm 2 x— 3x\nx)d x r+3r = 0, 

ix z In x \x z In x 
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cannot be algebrized systematically with this method, although it corresponds to 
pullbacks (not proper pullback) of differential equations with constant coefficients. 

In [21j . Bronstcin and Fredet developed and implemented an algorithm to 
solve differential equations over C(x,e^ f) without algebrizing it, see also [35]. As 
an application of proposition 12.3.21 we have the following resullQ. 

Proposition 2.3.5 (Linear differential equation over C(x, eSf), 0). Let f G C(z) 
be a rational function. Then, the differential equation 

% r -{f + !M_ fe f f a ^ ^ dxr + ^ ( e / /) ) 2 6 ( e / /) r = o, (2.6) 

is algebrizable by the change z = e$ t and its algebraic form is given by 

d 2 z y + a(z)d z y + b(z)y = 0, r(x) = y(z(x)). 
Proof. Assume that r(x) = y(z(x)), and z = z(x) = e^ dx . We can see that 

d x z = fz, d z y=^, ^=^(^)=^^(^-/+(^))^ 

replacing in d z y + a(z)d z y + b(z)y = we obtain equation (|2.6|) . □ 
Example. The differential equation 

figr - + i - 2xe x2 ^j d x r + \x 2 e x \ = 0, 

is algebrizable by the change z — e 2 and its algebraic form is given by 

d 2 y + 2zd z y + Xy = 0. 

Remark 2.3.6. According to proposition 1 2 . 3 . 51 we have the following cases. 

1. / = n^f^, for a rational function h, n G Z+, we have the trivial case, both 
equations are over the Riemann sphere and they have the same differential 
field, so that does not need to be algebrized. 

2. / = for a rational function h, n G Z + , (|2 .6[> is defined over an alge- 
braic extension of C(x) and so that this equation is not necessarily over the 
Riemann sphere. 

lr This result is given in [5] §2], but we include here the proof for completeness. 
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3- / =/= 9%^) f° r any rational function h, q 6 Q, (|2.6[) is defined over a transcen- 
dental extension of C(x) and so that this equation is not over the Ricmann 
sphere. 

To algebrize second order linear differential equations is easier when the term 
in d x r is absent, that is, in the form of equation (|1.2| and the change of variable 
is Hamiltonian. 

Definition 2.3.7 (Hamiltonian change of variable, [5]). A change of variable z = 
z{x) is called Hamiltonian if (z(x), d x z(x)) is a solution curve of the autonomous 
classical one degree of freedom Hamiltonian system 

d x z = d w H , w 2 

d x w = -d z H Wlth H = H(z,w) = -+V(z), 

for some V 6 K . 

Remark 2.3.8. Assume that we algebrize equation (|2.5[) through a Hamiltonian 
change of variables z = z(x), i.e., V € C(z). Then, K = C(z,d x z, . . .), but, we 
have the algebraic relation, 

{d x zf = 2h - 2V(z), h = H(z, d x z) € C, 

so that Ko = C(z, d x z) is an algebraic extension of C(z). By proposition 12.3.21 the 
identity connected component of the differential Galois group is conserved. On 
the other hand, we can identify a Hamiltonian change of variable z = z(x) when 
there exists a G K such that (d x z) 2 = a(z). Thus, we introduce the Hamilto- 
nian algebrization, which corresponds to the algebrization process done through a 
Hamiltonian change of variable. 

The following result, which also can be found in [51 §2], is an example of 
Hamiltonian algcbrization and correspond to the case of reduced second order 
linear differential equations. 

Proposition 2.3.9 (Hamiltonian Algcbrization, [5j). The differential equation 

d x r = q(x)r 

is algebrizable through a Hamiltonian change of variable z = z(x) if and only if 
there exist /, a such that 

— , - e C(z), where f(z(x)) = q(x), a(z) = 2(H - V{z)) = {d x z) 2 . 
a a 

Furthermore, the algebraic form of the equation d 2 r = q{x)r is 

dlv + \—d z y-^y = Q, r{x) = y(z(x)). (2.7) 
2 a a 
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Remark 2.3.10 (Using the Algebrization Method). The goal is to algebrize the 
differential equation d 2 r = q(x)r, so that we propose the following steps. 

Step 1 Find a Hamiltonian change of variable z = z(x) and two functions / and 
a such that q(x) = f(z(x)) and (d x z(x)) 2 = a(z{x)). 

Step 2 Verify whether or not f(z)/a(z) £ C(z) and d z a(z)/a(z) £ C(z) to see if 
the equation d 2 r = q(x)r is algebrizable. 

Step 3 If the equation d 2 r = q{x)r is algebrizable, its algebrization is 

<%V+ i— d z y - -y = 0, y(z(x)) = r(x). 
la a 

When we have algebrized the differential equation d 2 r = q(x)r, we study its 
intcgrability, Eigcnring and its differential Galois group. 

Examples. Consider the following examples. 

• Given the differential equation d 2 r = f(ta,nx)r with / £ C(tanx), we can 
choose z — z(x) — tanx to obtain ct(z) = (1 + z 2 ) 2 , so that z = z(x) 
is a Hamiltonian change of variable. We can see that ^«,^ £ <C(z) and 
the algebraic form of the differential equation d 2 r = /(tana;)r with this 
Hamiltonian change of variable is 

2z f(z) 

d2 * v + YTz^ dzV ~ (i + z 2 ) 2 y " 0, y( t&nx ) = r ( x )- 

• Given the differential equation 



dlr = 



2 



1 



we can choose z = z(x) = yl + x 2 to obtain 

f{z) = 2 ' a ( Z ) = 

z z z 

so that z = z{x) is a Hamiltonian change of variable. We can see that — £ 
<C(z) and the algebraic form for this case is 

1 „ z 2 + z-l 



d zy + -n — T\®*v 2 — 7— y = °: yWi + x 2 ) = r(x). 

Z(Z Z — I) 2^ — 1 

We remark that in general the method of Hamiltonian algebrization is not an 
algorithm, because the problem is to obtain a suitable Hamiltonian H satisfying 
definition 12.3.71 We present now a particular case of Hamiltonian algebrization 
considered as an algorithm^]. 



2 Proposition [2~3. 1 1 1 is a slight improvement of a similar result given in [5] §2], Furthermore, 
we include the proof here for completeness. 
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Proposition 2.3.11 (Hamiltonian Algebrization Algorithm, [5]). Consider q(x) = 
g{z\, ■ ■ ■ ,z n ), where Z{ = e XiX , A.j G C*. The equation d^r = q(x)r is algebrizable 
if and only if. 

l<i<n,l<j <n, geC(z). 

Furthermore, \ = c^A, where A € C* and Cj £ Q* and for the Hamiltonian change 
of variable 

\x p^ 

z = e~ , where c { = pi, qi e Z* , gcd(p,, = 1 and q = lcm(gi, ■ • • , q n ), 
the algebrization of the differential equation Ogr = q(x)r is 

o2 i 1q 29( zmi |"-i zm ") n IP* I I \\ I \ 

" z y+-o z y-q y = 0, rm = — , y\z[x)) = r{x). 

z X^-z 2 - qi 

Proof. Assuming Aj/Aj = cy £ Q* we can see that there exists A € C* and c, G Q* 
such that Ai = Ac;, so that 

e \ iX = e a\x = ei fAz = /A x \ q ' 5 Pu q t e z* 5 g C d(pi,^) = 1, lcm(<ji, ...,q n )=q. 
Now, setting 2 = z(a:) = e" 1 we can see that 

f(z) = g{z \...,z m, = 



Due to q|(7i, we have that m, G Z, so that 

d z a f 



G <C(z) 



and the algebraic form is given by 



9 z y + ~d z y - q 'r^-j y = 0, y(z(x)) = r(x). 

z A z z z 



□ 



Remark 2.3.12. Propositions 12.3.91 and 12.3.111 allow the algebrization of a large 
number of second order differential equations, see for example [77]. In particu- 
lar, under the assumptions of proposition 12.3. Ill we can algebrize automatically 
differential equations with trigonometrical or hyperbolic coefficients. 

Examples. Consider the following examples. 
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Given the differential equation 

fl2 e* a + Ze~i x -2ei x 1 2 5 A 1 x 3 

9 x r = eX + e - ix r > A i = 2' A2 = "3' As = 4' A4 = 1; As = "2' 

we see that S Q, A = 1, q = lcm(l, 2, 3, 4) = 12 and the Hamiltonian 

change of variable for this case is z — z{x) = eT5 x . We can see that 



\_ 2 . . _ z 6 + 3z- s -2z 15 d z a f 
144 z ' I( - z > ~ z i2 + z -is ' a ' a 



and the algebraic form is given by 

d 2 z y + -d z y - U4 Z + ^ ZiT V = 0, l/(en*) = r (s). 



• Given the differential equation 

d 2 x r = (e 2V2x + e-^ x - e 3x )r, Ai = 2^2, A 2 = -V2, A 3 = 3, 

we see that Ai/A 2 E Q, but A1/A3 ^ Q, so that this differential equation 
cannot be algebrized. 

We remark that it is possible to use the algebrization method to transform 
differential equations, although either the starting equation has rational coeffi- 
cients or the transformed equation has not rational coefficients. 

Examples. As illustration we present the following examples. 

• Consider the following differential equation 

2 x + 3x 2 — 5 

d ? r = 2T1 r = °> 

x z + 1 

we can choose z = z(x) = x 2 so that a = 4z and the new differential equation 
is 

<U+a^.y- 4z(z + 1) y = °. rt* ) = »■(*)• 

• Consider the Mathieu's differential equation 9 2 r = (a + bcos(x))r, we can 
choose z(x) = ln(cos(x)) so that a = e~ 2z — 1 and the new differential 
equation is 

d *y ~ i - e 2z d *y ~ °-\ y = °' 2/(M c °s(x))) = r ( x )- 

Recently, the Hamiltonian algebrization (propositions 12.3.91 and 12. 3. lip has 
been applied in [2j 02 [5] to obtain non-integrability in the framework of Morales- 
Ramis theory [68l [67] . 
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2.3.2 The Operator d z and the Hamiltonian Algebrization 

The generalization of proposition 12 .3 . 2| to higher order linear differential equations 
is difficult. But, it is possible to obtain generalizations of proposition 12.3.91 by 
means of Hamiltonian change of variable. We recall that z = z(x) is a Hamiltonian 
change of variable if there exists a such that (d x z) 2 = a(z). More specifically, if 
z = z(x) is a Hamiltonian change of variable, we can write d x z = ^fa, which leads 
us to the following notation: d z = \fad z . 

We can see that d z is a derivation because satisfy d z (f + g) = d z f + d z g and 
the Leibnitz rules 

d z (f- 9 ) = dj. 9+f .d z9 , djiy d *f-*-f-^ 



9/ g 2 

We can notice that the chain rule is given by d z (fog) = d g fogd z (g) ^ d g fogd z (g). 
The iteration of d z is given by 

d° z = i, d z = ^d z , d2 = V^d z d?- 1 = V^d z (...(V^d z )). 

•> v ' 

n times y/ad s 

We say that a Hamiltonian change of variable is rational when the potential 
V G <C(x) and for instance a £ C(x). Along the rest of this thesis, we understand 
d z = \fa.d z where z = z{x) is a Hamiltonian change of variable and d x z = ^/a. In 
particular, d z — d z = d x if and only if y/a = 1, i.e., z = x. 

Theorem 2.3.13. Consider the systems of linear differential equations [A] and [A] 
given respectively by 

d x Y = -AY, d z Y = —AY, A = [a l3 ], l=[a^], Y = [ Vil ], Y = [yn], 

where e K = <C(z(x),d x (z(x))), G C K = <C(z,y/a), 1 < i < n, 

1 < j < n j a ij( x ) = 0'ij{ z { x )) an d yn( x ) = yn( z ( x ))- Suppose that L and L are 
the Picard- Vessiot extensions of [A] and [A] respectively. If the transformation ip 
is given by 

x i — ► z 

^ ' yn{x)^> yn(z( x )) ' 

d x h-> d z 

then the following statements hold. 

• K~K, {K,d x )~(K,8 z )- 

• DG&\(L / K) ~ DGal(L / K) C DGal(i/C(z)). 

• (DGal(L/iT))° ~ (DGal(X/C(z))°. 
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• £{A)~£{A). 

Proof. We proceed as in the proof of proposition 12.3.31 As z = z{x) is a rational 
Hamiltonian change of variable, the transformation tp leads us to 

C(z)~cp(C(z))^K, K~K, C(z)^K, (A, d x ) ~ (A, d z ) 

that is, we identify C(z) with (p(C(z)), and so that we can view C(z) as a subfield 
of A and then, by the Kaplansky's diagram (see [531 1106] ). 




L K 




K 



C(z) 

so that we have DGal(L/A) ~ DGal(A/A) C DGal(A/C(z)), (DGal(A/A))° ~ 
(DGal(L/C(z)))°, and £{A) ~ £ (A). □ 

We remark that the transformation t^s, given in theorem 12.3.131 is virtually 
strong isogaloisian when y/a £ C(z) and for ^fa <E C(z), <p is strong isogaloisian. 
Furthermore, by cyclic vector method (see |81j). we can write the systems [A] and 
[A] in terms of the differential equations C and L. Thus, L is the proper pullback 
of C and£(£) -£(£). 



Example. Consider the system 



Oxl3 — e x +e -x 7l ~ e x +e -x 72) 

which through the Hamiltonian change of variable z = e x , and for instance \fa = z, 
it is transformed in the system 

&71 =--?$73. 

[A] := dJh = ^TTy73, 

2V2 ^ _ 2 2 -l 



9 Z 73 = 75+T71 - x (^ + i) 72- 
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One solution of the system [A] is given by 



1 



Z 2 + l 



z 
-z 



and for instance, 



1 




e 2x + 1 

is the corresponding solution for the system [A], 

Remark 2.3.14. The algcbrization given in proposition 12.3.91 is an example of how 
the introduction of the new derivative d z simplifies the proofs and computations. 
Such proposition is naturally extended to d x y + ad x y + by = 0, using ip to obtain 
&ly + cid z y + by = 0, which is equivalent to 

+ + ^) d z y + by = 0, (2.8) 

where y(x) — y(z(x)), a(z(x)) = a(x) and b(z(x)) = b(x). 

In general, for y(x) = y(z(x)), the equation F(d™y, . . . , y, x) = with co- 
efficients given by di k (x) is transformed in the equation F(d"y, . . . ,y, z) = 
with coefficients given by a !( .(z), where ai k (x) = ai k (z(x)). In particular, for 
\fa, ai k £ C(z), the equation F(d™y, . . . ,y,z) = is the Hamiltonian algebrization 
of F(d%y, . . . ,y,x) =0. Now, if each derivation d x has order even, then a and a,i k 
can be rational functions to algebrize the equation F(d™y, ...,y,x) = 0, where 
ai h £ <C(z(x), d x z{x)). for example, that happens for linear differential equations 
given by 

d 2 x 1l y + a n - 1 {x)dl n - 2 y + ... + a 2 (x)d^y + ai{x)d 2 x y + a (x)y = 0. 

Finally, the algebrization algorithm given in proposition 12.3.111 can be naturally 
extended to any differential equation 

F(d"y,d2- 1 y,...,d x y,y,e» t )=0, 

that by means of the change of variable z = e^ ix is transformed into 

F(^y,^- 1 y,...,9 z y, 2 /,z) = 0. 

Particularly, we consider the algebrization of Riccati equations, higher order linear 
differential equations and systems. 
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Examples. The following corresponds to some examples of algebrizations for dif- 
ferential equations given in [HH p. 258, 266]. 

1. The equation C := d 2 y + (—2e x — l)d x y + e 2x y = with the Hamiltonian 
change of variable z = e x , \fa = z,a = — 2z — 1 and b = z 2 is transformed in 
the equation £ := d 2 y — 2d z y + y = which can be easily solved. A basis of 
solutions for C and L are given by (e z , ze z ) and (e e , e x e e ) respectively. Fur- 
thermore K = C(e x ), K = C(z), L and L are the Picard-Vessiot extensions 
of C and £ respectively. Thus, DGai(L/K) = DGal(L/A). 

2. The differential equation 

. „ 2 -24e*-25„ 20e* 

£ := o^w H a x v H y = 

xU Ae x + 5 y 46^ + 5 

with the Hamiltonian change of variable z = e x , *J~u = z, 

^ -24z- 25 - 20z 
and b 



4z + 5 4z + 5 

is transformed in the equation 

£ := ^ + x(4z + 5) d ' y + ^4zT5) y = °' 

which can be solved with Kovacic algorithm. A basis of solutions for £ is 
(z + l,z 5 ), so that a basis for £ is (e x + l,e 5x ). Furthermore K = C(e x ), 
K = C(z), L and L are the Picard-Vessiot extensions of C and C respectively. 
Thus, DGal(L/A') = DGal(£/A) = e. 

Remark 2.3.15 (Algebrization of the Riccati equation). The Riccati equation 

d x v = a(x) + b(x)v + c(x)v 2 (2.9) 

through the Hamiltonian change of variable z — z(x), becomes in the Riccati 
equation 

d z v = -^=(a{z) + b{z)d + d(z)v 2 ), (2.10) 

where v(x) = v(z(x)), a{z{x)) = a(x), b{z{x)) = b(x), c(z{x)) = c(x) and 
y/a(z(x)) = d x z(x). Furthermore, if yfa, a, b, c € C(x), equation (|2.10| is the 
algebrization of equation (|2.9p . 

Example. Consider the Riccati differential equation 

C := d x v = [ tanh a; ) v + (3 tanh a; — 3 tanh 3 x) v 2 , 

\ tanh x ) 
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which through the Hamiltonian change of variable z = tanhx, for instance y/a ■ 
1 — z 2 , is transformed into the Riccati differential equation 

C := d z v = — v + 3zv 2 . 

z 



One solution for the equation £ is 

1 



, being c a constant , 



z(3z — c) ' 

so that the corresponding solution for equation C is 

1 



tanh x(3 tanh x — c) 

The following result is the algebrized version of the relationship between the 
Eigenrings of systems and operators. 

Proposition 2.3.16. Consider the differential fields K , K and consider the systems 
[A] and [A] given by 

d x X = -AX, d z X = -AX, d z = V^d z , A = [ay], A = [ciij], a tJ G K, 3y G K, 

where z = z(x), X(x) = X(z(x)) 7 ciij(z(x)) — aij(x), then £(A) ~ £(A). In 
particular, if we consider the linear differential equations 

n—l n—1 

C:=d:y + Y / a kd k x y = and Z := d n z y + £ a k d k z y = 0, 

k=0 k=0 

where z = z(x), y(x) = z((x)), a k {z{x)) = ak(x), a k G K, ctfc G K, then £ (£) ~ 
£(£), where C := Zy — and C := £/y = 0. Furthermore, assuming 



/ 1 ... 

Pnl ■■■ Pnn) \— ^0 _a l ••■ —^n-Xj 



Pll ■■■ Pin 
■ 



then 

£{Z) = J Y.Pikd'lT 1 : d x P = PA- AP, p ik G K 1 , 



.fc=i 



if and only if 

£(£) = J Y.P^' 1 : = PA-AP, p lk G A' 



,fc=i 
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Proof. By theorem EXH we have that K ~ K, 8(A) ~ £(1) and £(£) ~ 8(2). 
Using the derivation d z and by induction on lemma l2.3.3l we complete the proof. 

□ 

Examples. We consider two different examples to illustrate the previous proposi- 
tion. 

• Consider the differential equation L\ := d 2 y — (1 + cosx — cos 2 x)y = 0. 
By means of the Hamiltonian change of variable z = z(x) — cosx, with 
\fo, = — vl — z 2 , Ci is transformed into the differential equation 

? «2~ z a - 1 + 2- 

£i := c£y - ^d z y — y = 0. 

1 — z z 1 — z z 

Now, computing the eigenring of Z\ we have that £(£i) = Vect(l), therefore 
the eigenring of £i is given 8 (£i) = Vect(l). 

• Now we consider the differential equation £2 := d 2 y = (e 2x + |j y. By means 
of the Hamiltonian change of variable z = e x , with s/a = x, £2 is transformed 
into the differential equation 

C 2 :=d 2 y+- z d z y-(l + ^y = 0. 
Now, computing the eigenring of £2 wc have that 

S(2 2 ) = Vect (l, -2 (f^I) 9, - = Vect (l, -2 (f^I) & - ^) , 

therefore the eigenring of £2 is given by 

W = Vcct (l. - 2 (^) * - ^) ■ 
The same result is obtained via matrix formalism, where 

A=(j + ^ ^,3=^° 9 j),d x P = PA-AP,aP = PA-AP, 
with P e Mcrf(2,C(e x )) and P e Mai(2,q»). 

2.3.3 Applications in Supersymmetric Quantum Mechanics 

In this subsection we apply the derivation d z to the Schrodingcr equation H^f = 
A^P, where H = —d 2 + V(x), V £ K . Assume that z = z(x) is a rational Hamil- 
tonian change of variable for = X^>, then K = C(z(x), d x z(x)). Thus, the 
algebrized Schrddinger equation is written as 

ff$ = A$, H = -d 2 + V(z), d 2 =adl + ^d z ad z , K = C(z,Va). (2.11) 
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The reduced algebrized Schrodinger equation, obtained through equation (|1.3p . is 
given by 

V(*)=V + ^, 

H$ = A$, H = (-<9 2 2 + V(z)) , V = d z W + W 2 , (2.12) 

v\; — I Wf) 

^ — 4 a(z) ' 

The eigenfunctions ^P. and <f> corresponding to the operators H, H and H arc 
related respectively as 

$(z(x)) = yofif(z(x)) = ^a*(a;). 

In order to apply the Kovacic's algorithm we only consider the algebrized 
operator H, whilst the eigenrings will be computed on H. Also it is possible to 
apply the version of Kovacic's algorithm given in reference |101j to the algebraized 
operator H . The following results are obtained by applying Kovacic's algorithm 
to the reduced algebrized Schrodinger equation (equation (|2.12[) ) H$ = A4>. 

Proposition 2.3.17. LetL\ be the Picard-Vessiot extension of the reduced algebrized 
Schrodinger equation H$ = A<!> with a, V G C[z]. If deg a < 2 + degV , then 
DGal(L/A') is a not finite primitive for every A G A. 

Proof. Suppose that deg a = n and deg V = m. The reduced algebrized Schrodinger 
equation H$ = A$ can be written in the form 

~ ^ Aad 2 a — 3(d z a) 2 + 16a(V - A) 
d 2 z <S> = r$, r = z - / 1 '-. 

Due to m > n — 2 we have that o(r oc ) = n — m < 2, which does not satisfy the 
condition (oo) of the case 3 of Kovacic's algorithm, therefore DGal(L/A") is a not 
finite primitive for every A 6 A. 

□ 

Proposition 2.3.18. Let La be the Picard-Vessiot extension of the reduced algebrized 
Schrodinger equation H$ = \<S> with a G C[z], V € C(z). If o(V)oo < 2 — deg a, 
then DGal(L,\/A") is a not finite primitive for every A G A. 

Proof. Suppose that V = s/t, being s and t co-primes polynomials in C(z). Assume 
that deg a = n, degs = m and degt = p. The reduced algebrized Schrodinger 
equation H$ = A$ can be written in the form 

4tad 2 a - 3t(d z a) 2 + I6a(s - \t) 

dt§ = r$, r= - — -. 

16ta 2 
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Due to m > n+p — 2 we have that o(r oc ) = p + n — m < 2, which does not satisfy 
the condition (oo) of the case 3 of Kovacic's algorithm, therefore DGal(L/A") is a 
not finite primitive for every A e A. 

□ 



Remark 2.3.19. In a natural way, we obtain the algebrized versions of Darboux 
transformation, i.e., the algebrized Darboux transformation, denoted by DT. By 
DT n we denote the n iteration of DT, and by CI n we denote the algebrized Crum 
iteration, where the algebrized wronskian is given by 



2/1 



W(y±, ...,y„) 



ill 



Vn 



bn-l~ 



In the same way, we define algebrized shape invariant potentials, algebrized super- 
potential W, algebrized super symmetric Hamiltonians H±, algebrized supersym- 

^ ^ ( — ) — f -^—dz 

metric partner potentials V±, algebrized ground state VPq — e J 73 , algebrized 
wave functions \E'^ , algebrized raising and lowering operators A and A^ . Thus, 
we can rewrite entirely the section 12.21 using the derivation d z . 



The following theorem show us the relationship between the algebrization 
and Darboux transformation. 



Theorem 2.3.20. Given the Schrddinger equation £\ := H_^/^ ' = >, the 
relationship between the algebrization tp and Darboux transformations DT, DT 
with respect to C\ is given by DT</? = ^DT, that is DTy>(£) = ipDT(C). In 
other words, the Darboux transformations DT and DT are intertwined by the 
algebrization tp. 

Proof. Assume the equations^ := = A* (_) , Z\ := iL§ ( ~) = 

£\ := = A*<+) and Z := = A$M, where the Darboux transfor- 

mations DT and DT are given by DT(£) = Z, DT(£) = Z, 

V- V-^ V+ 

DT : , DT : 

and tp(C\) = £\, where the algebrization tp is given as in theorem 12.3.131 Then 
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the following diagram commutes 



r DT - 7 
L\ *- C 



DT 



BT(p(£) = ipDT(C) & C = C. 



□ 



To illustrate DT we present the following examples. 



Examples. Consider the algebrized Schrodinger equation > = A^ with: 



yja{z) = Vz 2 - 1 and V-(z) = j^. Taking Ai = 1 and W\ 
have that DT(V_) = V+(z) 



(-) 



z+l 



and 



we 



DT($^ ') = * 



;(+) 



= ^ 



1 



(-) 



where ^ is the general solution of H-ty^ ' = A^( ) for A ^ 1. 

The original potential corresponding to this example is given by V_ (x) = 
c ° sh - 1 and for Ai = 1 the particular solution is given by sl , nh - . . 

cosnx- 1 1 o ^ cosh a:— 1 

Applying DT we have that DT(V_) = V+(a;) = ggf^ and DT(*^ ) ) = 



( + ) 

A 



(-) 



• V a = 

DT(V_ 



A ^ sinh x A 
2 



(-) 



, V_(z) = z 2 - z. Taking Ai = and \I> 
: V+ = z 2 + z and DTY*!"') = $ 



(-) 

A 



( + ) 



(-) 




e z we have that 
z^^\ where 



is the general solution of H_^>( 



corresponds to the Morse potential V_(x) 
list (TOPI) , 



A ^ z A 

A$(-) for A ^ 0. This example 
= e~ 2:c — e~ x , introduced in the 



To illustrate CI„ we present the following example, which is related with the 
Chebyshev polynomials. 



Example. Now, considering \fot = —yl — z 2 , V — with eigenvalues and eigen- 
functions Ai = 1, A 2 = 4, $i = z, $ 4 = 2z 2 - 1, $„2 = T n (z), where T„(z) is 
the Chebyshev polynomial of first kind of degree n. The algebrized Wronskian for 

n = 2 is 

W(z, 2z 2 - 1) = -\J\ - z 2 (2z 2 + 1), 
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and by algebrized Crum iteration we obtain the potential 

Gl 2 {V) = V 2 = {{2z 2 - l)d 2 z + zd z ) In VF(z, 2z 2 - 1) 
and the algebrized wave functions 

- (2) _ W(z,2z 2 -l,T n ) 



CI 2 (*a) = *\ 



W(z,2z 2 - 1) 



In a natural way we introduce the concept of algebrized shape invariant po- 
tentials V n -\-i_{z, a n ) — V n (z, a n +i) + i?(a„), where the energy levels for n > are 
given by E n — R(ai) + ■ ■ ■ R(a n ) and the algebrized eigenfunctions are given by 
^n(di) = A'(z, ai) ■ ■ ■ A'(z, a n )^o{z, o, n ). To illustrate the algebrized shape in- 
variant potentials and the operators A and A\ we present the following example. 

Example. Assume ^fa. = 1 — z 2 and the algebrized super potential W(z) = z. 
Following the method proposed in remark [2.2.161 step 1, we introduce /i G C to 
obtain W(z; /i) = fiz, and 

V ± (z;fi)=W 2 (z;fi)±d z W{z;Li)=^Tl)z 2 ±^ 

thus, V+(z; ao) = ao(ao — l)z 2 + ao and V-(z; ai) = ai(a\ + \)z 2 — a\. By step 2, 

d z (V + (z;a ) - VL(«;oi)) = 2z(l - z 2 )(a (a - 1) - 01(01 + 1)). 

By step 3, we obtain 

ai(ai + 1) = a {a - 1), a\ - a\ = -(at + a ), 

and assuming a\ ^ ±ao we have a\ = /(ao) = ao — 1 and R(ai) = 2ao + 1 = 
(ao + l) 2 — ap = a 2 — a 2 . This means that the potentials V± are algebrized shape 
invariant potentials where E — E n is easily obtained, 

n n 

E n = ^2 R ( a k) = J! ( a k ~ a k-l) = a n - a 0> «» = /"(°o) = a + Tl. 
k=l fc=l 

Now, the algebrized ground state wave function of V— {z, ao) is 

$ = e J 



Finally, we can obtain the rest of eigenfunctions using the algebrized raising op- 
erator: 

ao) = A\z, oq)A\z, a\) ■ ■ ■ jv(z, a„_i)* ( 
This example corresponds to Poschl- Teller potential introduced in the list l| 1.401) . 
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Now to illustrate the power of Kovacic's algorithm with the derivation d z , 
we study some Schrodingcr equations for non-rational shape invariant potentials 
given in list (|1.40|) . We work with specific values of these potentials, although we 
can apply our machinery (algcbrization method and Kovacic's algorithm) using 
all the parameters of such potentials. 

Morse potential: V(x) = e~ 2x — e~ x . 

The Schrodinger equation = A^P is 

fig* = (e~ 2x - er x - A) #. 

By the Hamiltonian change of variable z = z(x) = e~ x , we obtain 

Z 2 - Z - i 

a(z) = z 2 , V{z) = z 2 -z, V(z) = i. 

Thus, K — C(z) and K = C(e x ). In this way, the algebrized Schrodingcr equation 
H $ = A$ is 

z 2 d 2 $ + z<9 2 $ - (z 2 - z - A)$ = 
and the reduced algebrized Schrodinger equation H$ = A<£> is 

o z 2 z -7 A 



di<& = r$, r 

z~ 

This equation only could fall in case 1, in case 2 or in case 4 (of Kovacic's algo- 
rithm). We start analyzing the case 1: by conditions and 003 we have that 



[Vr\ = 0, ao= 1±2 2 A , [V?L = 1. and Q t,=^\- 

By step 2 we have the following possibilities for n G Z + and for A £ A: 

,2 



A++) 


/V 


= a to 


-a+ = 


-1-V 


-A, 


A = 


— (n 


A+-) 


n 


= a to 


-olq = 


-l + V 


=A, 


A = 


— (n 


A-+) 


n 


= a ^o 


- "0 = 






A = 


—n, 


A—) 


n 


= a ^o 


- "0 = 






A = 


—n. 



We can see that A £ A_ = {— n 2 : n E Z + }. Now, for A 6 A, the rational function 



84 



Chapter 2. Differential Galois Theory Approach to SUSY QM 



uj is given by: 



A++) 


LU = 


1 


+ 


3+2n 
2z ' 


A e A++, 


r n = 


4n 2 +8n+3 , 
4z 2 + 


2n+3 , 

z 


1, 


A+-) 


LU = 


1 




l+2n 

2z ' 


A e A+_, 


r n = 


4n 2 +8n+3 
4z 2 


2n+l 

z 


1. 


A-+) 


LU = 




1 


i l+2n 
2z ' 


A 6 A_+, 


r n = 


4n 2 -l 2n-\ 
4z 2 z 






A—) 


LU = 




1 


, l-2n 
+ 2z ' 


A 6 A__, 


r n = 


4n 2 -l , 2n- 
4z 2 1 z 







where r n is the coefficient of the differential equation c^<f> = r„<I>. 

By step 3, there exists a polynomial of degree n satisfying the relation (|1.6|) . 
A ++ ) 9, 2 P„ + 2(l + 2±f»)^P n + ^Pn = 0, 
A+_) dlP n + 2 (1 - ±±fa) 5 Z P„ + H^l p„ = o, 

A_+) 9fP n + 2(-l + i±f^)a z P n + ^ i P„ = 0, 
A__) 5?P„ + 2 (-1 + i=fa) d z ?n + ^P n = 0. 

These polynomials only exists for n = A = 0, with A € A |_ U A So that 

the solutions of H*f> = 0, = and H$ = are given by 

$o = \/ze~ z , § = e~ 2 , * = e~ e . 

The wave function 'J'o satisfy the conditions (|1.23[) . which means that is ground 
state (see [31]) and € spec p (iJ). Furthermore, we have 

DG&l(L /K) = BGal(L /K) = DGal(L /C(z)) = B, 

£{H) = £{H) = £(H) = Vect(l). 

We follow with the case two. The conditions C2 and 003 are satisfied, in this 
way we have 

E c = |2,2 + 4\/^A,2-4V^A| and E oo = {0}, 

and by step two, we have that 2 ± \/— A = m G Z + . so that A = — (j 1 ^-) 2 and 
the rational function 6 has the following possibilities 



z 
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By step three, there exist a monic polynomial of degree m satisfying the recurrence 
relation p.7[) : 

Q^P -f- 3m+6 g2 jj 4z 2 — 4z — 2m 2 — 7m— 6 g 4mz+8z- 4m— 6 jj _ q 

e-) d 3 z P m - ^fd 2 z p m - **'- 4 *-? m '- m d s p m + Amz -jr i - 2 ? m = o. 

We can see that for m = 1 the polynomial exists only for the case 6L, being 
Pi = z — 1/2. In general, these polynomials could exist only for the case 9- with 
m = In — 1, n > 1, that is A £ {— rt 2 : n > 1}. 

For instance, by case one and case two, we obtain A = {— n 2 : n > 0} = 
speCp(P). Now, the rational function cj> and the quadratic expression for u) are 



In - 1 d z P 2n -i 9 -M ± VM 2 - 47V 
= + z - 2n , w 2 + Mw + A r = 0. w = , 

z P 2 n-1 2 

where the coefficients M and TV are given by 

M _ 2n-l gA-l N _ n 2 -n+\ ( 2n ~ ^^ife 1 ~ 2 | c\ 2 P 2n -i 
z Pln-l 2,2 z Pln-\ 



Now, A = M 2 - 47V ^ 0, which means that H$ = -n 2 <3> with n <G Z+ has two 
solutions given by Kovacic's algorithm: 

$l,n = " > $2,n .= " ■ 



The solutions of = — to \t are given by 



o _P n e~ z $ P n _ ie ' 

Wl„ - " , W 2 ,n 



Z" Z" 

and therefore, the solutions of the Schrodinger equation H^f = —n 2 ^ are 

*i,n = P n e- e ~°e nx , * 2 .„ = P„_ie e ~ V 1 *, P n = P n oz. 

The wave functions ^\. n = satisfies the conditions of bound state, and 
for n — 0, this solution coincides with the ground state presented above. Therefore 
we have 

$n = $0fnPn, *n = *o/A f n (z) = — . 

Z 

Thus, the bound states wave functions are obtained as 

*„ = *o/„P«, U(x) = f n (e- x ) = e nx . 
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The Eigcnrings and differential Galois groups for n > satisfies 

DGal(L n /K) = DGal(L n /K) = DGal(L„/C(z)) = G m , 

dim c £(H + n 2 ) = dim c £(H + n 2 ) = dim c £(H + n 2 ) = 2. 

We remark that the Schrodinger equation with Morse potential, under suit- 
able changes of variables [60], falls in a Bessel's differential equation. Thus we can 
obtain its intcgrability by means of corollary II. 1.321 

It is known that Eckart, Rosen-Morse, Scarf and Poschl- Teller potentials, 
under suitable transformations, fall in an Hypergeomctric equation which allows 
apply theorem 11.1.301 These potentials are inter-related by point canonical co- 
ordinate transformations (see [24J p. 314] ), so that A = C due to Poschl- Teller 
potential is obtained by means of Darboux transformations of V = ( 66, 86 ). 
We consider some particular cases of Eckart, Scarf and Poschl- Teller potentials 
applying only the case 1 of Kovacic's algorithm. The case 1 allow us to obtain the 
enumerable set A„ C A, which include the classical results obtained by means of 
supersymmetric quantum mechanics. Cases 2 and 3 of Kovacic algorithm also can 
be applied, but are not considered here. 

Eckart potential: V(x) = 4coth(x) + 5, x > 0. 

The Schrodinger equation = X^> is 

d 2 x -$> = (4coth(x) + 5- A)*. 

By the Hamiltonian change of variable z = z(x) = coth(a;), we obtain 

a(z) = (l-z 2 ) 2 , V>(z) = 4z + 5, V(z)= 

Thus, K — C(z) and K = C(coth(rr)). In this way, the algebrized Schrodinger 
equation H'fy = A\& is 

(1 - z 2 ) 2 d 2 $ - 2z(l - z 2 )<9 z $ - (4z + 5 - A)$ = 

and the reduced algebrized Schrodinger equation H$ = A$ is 

9 2 $ = r$, 



4z + 4-A 2-A |_l _A l_| 



{z-l) 2 {z + lf {z-lf {z-1) ( z + i) 2 (z + l) 

We can see that this equation could fall in any case of Kovacic's algorithm. 
Considering A = 0, the conditions {ci,C2,ooi} of case 1 arc satisfied, obtaining 

Wr]-i = [Vr]i = [Vr]oo = a+ = 0, at t = a~ = 1, = 2, a~{ = -1. 
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By step 2, the elements of D are and 1. The rational function ui for n = and 
for n = 1 must be 

1 -1 

LU = 1 . 

z + 1 z - 1 

By step 3 we search the monic polynomial of degree n satisfying the relation (|1.6p . 
Starting with n = the only one possibility is Pq(z) = 1, which effectively satisfy 
the relation (|1.6|) . while -Pi(z) = z + ao does not exists. In this way we have 
obtained one solution using Kovacic algorithm: 



2+1 - 2+1 
$o = 7, = 



2-1' u y(z-l) 3 ' 
this means that S A„. We can obtain the second solution using the first solution: 

z 2 + z - 4 - 4 In (2 + 1) 2 - 4 In (2 + 1) - $0,2 
$0,2 = : , *0,2 = 



2-1 ' Vi^T 

Furthermore the differential Galois groups and Eigcnrings for A = are 
DGal(L /C(z)) = G OJ DGal(L /K) = DGal(£ /X) = G {2} , 

£(H)=Vect(l,£^ ■ 2(2 + 1) 



(1-z) 2 - {l-zf 

= VCCt I, 1 ' (1^)2^ - (1 - ,)3 

= Vect (l, (coth(,) + l) 2 _ coth 2 (,) + 3coth(,) + 2 

V (l-coth(a;)) 2 (l-coth 2 (a;)) (1 - coth(x)) 3 

Now, for A 7^ 0, the conditions {02, 001} of case 1 arc satisfied: 
Wr\-i = [\/r\i = [\fr\oo = a+ = 0, a x = 1, 

— 1 2 ' 1 2 

By step 2 we have the following possibilities for n G and for A G A: 

A++-) n = a+ - a\ -a^ = -l- vT=*->/9=* A = 4 - - n 2 - 2n, 

A+— ) n = a+ - all - af = -1 + v^+v^ , A = 4 - ^pjp- - n 2 - 2n, 
A_+_) n = a~ -a^-ar= ^+^ , A = 5 - £ - n 2 , 

A—) n = a--aZ 1 -ai=>£EE±V9EE j A = 5 - 4, - n 2 . 
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Therefore, we have that 

A„ C j 4 - - — — ^ - n 2 - 2n : n £ Z + 1 U j 5 - 4r - n 2 : n £ Z + 
^ [n+ 1) J I n 

Now, for A G A, the rational function a; is given by: 

a \ _ z(n~l)-re 2 -2ra-l _ -2z 2 (n- l)+4z(ra+l) 2 + (n + l) (ra 3 +3n 2 + 2n + 2) 

A + + _J W — ( n + i)( z + i)(,_l) I r„ — (n + l^Oz+l^z-l)* ' 

* \ _ 7iz(n + l) + 2 _ rez 2 (n+l) 3 +4z(n + l) 2 + re 3 +2re 2 +n+4 

A + 1 U ~ (n + l)(z + l)(l-z) ' Tn ~ ( n + l)2( 2 + 1 )2( z _ 1) 2 , 

. \ _ z(n-2)-n 2 _ -2z 2 (n-2) + 4re 2 z + n(n 3 -?i+2) 

+ U — „( Z + 1 )( Z _ 1 ) , r n — n2(z + l)2( z _l)2 ) 

A \ _ nz(n-l) + 2 _ n 3 z 2 (ra-l) + 4re 2 z + n 3 -n 2 +4 

A J W — n(«+l)(l-») ' r " _ n 2 (z+l)2(z-l)2 ' 

where r„ is the coefheient of the differential equation d 2 Q = r n $. 
By step 3, there exists a monic polynomial of degree n satisfying the relation (|1.6p . 

A \ r) 2 P | o I z(n-\)—n 2 -2n-\ \ o p , 2(1- ra) p n 

^M-+-J w 2 f n T z ^ (n+l)(z+l)(z-l) J + ((ri+l) 2 (z+l)(z-l) ^™ ~~ U ' 

A 1 r) 2 P 4- 9 ( nz(n+l) + 2 \ o p , n(ra+l) p — fl 

J1 + / U z r ™ + Z \(n+l)(z+l)(l-z) J u z r n T (^ + l)( z _l) r n ~ u < 

A . 1 fl 2 P + 2 f ^(»- 2 )-" 2 ^ a B j 2 ( 2 ~») p _ n 

The only one case in which there exist the polynomial P n of degree n is for A^ ) . 

The solutions of the equation H$ = A<t>, with A ^ 0, are: 

A++-) $„=P n /„$0, $0 = ^ /„= (Z-1)*+If(z + 1)^+If, 

^ ^ , , n(l-n) -n(n + 3) 

A+__) $„ = P„/„$ , $0 = ^ /n = (z-l)WTT(z + l)^+TT, 

A_+_) $ n = P n / n $ 1 , $1 = ^ /„ = (z + l) s ^ 1 (z-l) "" 2 t 3 "" 2 , 

A—) $„ = P„7„*1, $1 = §_^ In = + l)" 2 ^ 2 (* ~ 1) ~" 22 """ 2 ■ 

In any case = y-r^ , but the case A_| ) includes the classical results obtained 

by means of supersymmetric quantum mechanics. Thus, replacing z by coth(:r) we 
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so 



obtain the eigenstates The Eigenrings and differential Galois groups for n > 
and A G A„ satisfies 

DGal(A A /A") C G {2m} , DGal^/AQ Q G {2m} , 

DGal(L A /C(z)) = G m , 
dim c(z) £(H + A) = 2, £ (# + A) = £{H + A) = Vect(l). 



Scarf potential: V(a;) 



sinh x — 3 s inh x 



cosh x 

The Schrodinger equation Hty = is 

„ 9 T / sinh 2 a; — 3 sinh a; A 

d 2 * = J E *. 

V cosh a: / 

By the Hamiltonian change of variable z = = sinh(a;), we obtain 
a(z) = l + z 2 , V(z) 



1 + z 2 



Thus, A" = C(z, y/l + z 2 ) and AT = C(sinh(a;), cosh(x)). In this way, the reduced 
algebrized Schrodinger equation H$ = A<J> is 

^( ^Z 1 ?^- 1 )** ^3-4A. 

2 V 4(z 2 + i) 2 y 

Applying Kovacic's algorithm for this equation with A = 0, we see that does 
not falls in case 1. We consider only A ^ 0. By conditions {02,002} of case 1 we 
have that 



[Vrl-i = [Vrli = [Vrloo = o, a ±=±±^, 

-1 4 2' — 1 i~ 2 1 4 ' 2> 1 4 2' 

By step 2 we have the following possibilities for n € Z + and for A G A: 
A +++ ) n = a+ - ati - = ^Y" 4 , A = 4?i 2 + 16n + 15, 

A+__) n = a+ - QQ - a t ~ = ^ y + T+2 , A = 4n 2 - 8n + 3, 
obtaining in this way 

A„ C {4n 2 + 16n+ 15 : n £ Z+} U {4n 2 - 8n + 3 : n £ Z+} . 
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Now, the rational function u> is given by: 

A +++) " = --5—— , A + __) 



2(z 2 + l) 2(z 2 + l) 

By step 3, there exists Pq = 1 and a polynomial of degree n > 1 should satisfy 
cither of the relation (|1.6[) . 

A "\ fl 2 P _1_ 5z-2 a p rtz 2 (n+4)+3z+« 2 +4n-3 p _ n 
A + + +J + -£T^lO z f n (« a +l) a ™ — U ' 

A ~\ fl 2 P i 2-2 a p nz 2 (Tt-2)+32+n 2 -2rt~3 g _ n 
-<H 1 °z r n T .i +1 t/;r n (p+ip " ~~ 

In both cases there exists the polynomial P n of degree n > 1. Basis of solutions 
{$i n , $2,n} of the reduced algebrized Schrodinger equation are: 

A+++) = Pn/n*l,0, *1,0 - (1 + z 2 ) § e~ arctan z , /„ = 1, 

$2,„ = <9n?n$2,0, $2,0 = 22 + 21 |+gW e ~ «°t~» ^ = L 

A+— ) *!,„ = P„/n$l,0, $1,0 = ^7=f e arctan *, /„ = 1, 

^ rh rh 1 „arctan z P 1 ^ — 2 arctan g ^„ C: i 

3 > 2,n — Vn5n*2,0> * 2 ,0 — 4/ 1+z ^ e J ^ 1+z2 e dz i 9n — 1. 

In both cases = ^JL a , but the classical case (see references [231 El) is 
A_| ), so that replacing z by sinhx and A by 3 — 4E we obtain the eigenstates 

The Eigenrings and differential Galois groups are 

£(H - A) = £(# - A) = £(H - A) = Vcct(l), 
DGal(L x /K) = DGal(£ A /X) = DGal(L A /C(x)) = B. 



Poschl- Teller potential: V(t) = ^^^i^^giJ^].^ 2 > x > 0. The reduced alge- 
brized Schrodinger equation H<E> = i?$ is 



: r,' "'V - )*, A = 3-4£. 



Az 4 - (A + 3)z 2 + 
4z 2 (z 2 - l) 2 

Considering A = and starting with the conditions {02,001} of case 1, we 
obtain 

[y/r\o = Wr\-i = iVr\i = [VrU = a+ = 0, a" = 1, 
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By step 2, the elements of D are and 1. The rational function to has the following 
possibilities for n = and for n = 1: 



A++-) n = 0, c = ^ + ^ + 

A_| | ) n = 0, uj 

A | ) n = 1, tj 

A | ) n = 1, 

By step 3 we search the monic polynomial of degree n satisfying the relation (|1.6|) . 

Starting with n = the only one possibility for A + _| ) and A_| | ) is Pq(z) = 1, 

which does not satisfy the relation (|1.6[) in both cases, while Pi(z) = z + ao 

effectively does exists, in where ao = — § for A | ) and ao = | for A | ). In 

this way we have obtained two solutions ($i,o, < &2,o) using Kovacic's algorithm: 



5/4 


-1/4 


Z+l T 


z-1 


-1/4 


5/4 


z+l 


r z-1 


5/4 


-1/4 


z+l 


z-1 


-1/4 


5/4 


z + l 


r z-1 



*i.o = (i-^) T-^ti *i.o = C 1 - c) ^r. 



*>.o = (l + ^)y^. * 2 .o = (l + 
this means that £ A n . Furthermore, 

DGal(L /C(x)) = G [4] , DGal(Z /A") = DGal(L /K) = e, 

dim c £(H) = 2, dim c £ (#") = dim c 5(F) = 4. 
Now, for A 7^ we see that conditions {c2,ooi} of case 1 leads us to 

[yffio = Wr\-! = [Vr]i = [VfU =0, a± = 
aj: 1 =af = |, al 1 = a 1 " = -i, = 2, = -1. 
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By step 2 we have the following possibilities for n £ Z + and for A 6 A: 



A++++ 
A+++- 
A++- + 
A++- 
A+-++ 
A+-+- 
A+— + 
A+— 



n = "oo - a -l ~ a l - a Q 



n = a' — a i — av — a n = 



n = a+ — a_j — a x — 



n = a ' -a , - a, — a n = 



n = Q!+ — a i — at — ott 



n = a oo- a -l - a l - a 6 = 



+ _ VA+T-8 



2 

/X+T-2 



2 

/ A+T-5 



2 



2 

^+T-5 



2 

/ A+T+1 



2 

^+T-2 



2 

/A+T+4 



obtaining A„ C A a U A h U A c U A d U A e , where 



A = An 2 + 32n + 63, 
A = 4n 2 + 8n + 3, 
A = 4n 2 + 20n + 24, 
A = 4n 2 - in, 
A = 4n 2 + 20n + 24, 
A = 4n 2 - in, 
A = 4n 2 + 8n + 3, 
A = 4n 2 — 16n+15, 



A a = {4n 2 + 32n + 63 : n G Z + } , A 6 = {4n 2 + 8n + 3 : n e Z + } 
A c = {4n 2 + 20n + 24 : n G Z+} , A d = {4n 2 - in : n G Z+} , 
A e = {4n 2 - 16n + 15 : n G Z + } . 



Now, the rational function lo is given by: 



+++ 



5/4 ,5/4,2 



> u - 7+1 + — 



A 



+++- 



A++- + J u= z 
A+-++) u> = 



_ 5/4 , -1/4 , 2 



-1 



— 



-1/4 , 5/4 

Z+l "T" Z-l 



A+— +) 



-§, A+_+_) w 

-144 ^V4 2 A+ ) W 



5/4 , 


5/4 , 


-1 


z + l T 


z-l 


z ' 


5/i 4 




-1 


z + l T 


z-l 1 


z ' 


-1/4 


f s/f. 


-1 


z+l 


r z-l n 


z ' 


-1/4 


f - 1/4 




z+l 


r z-l 





By step 3, there exists a monic polynomial of degree n satisfying the relation 
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(fl~6l) . 

A++++ 
A+++- 
A++-+ 
A++— 
A+-++ 
A+-+- 
A+— + 
A+ 

The polynomial P n of degree n exists for A„ 6 A5 with n even, that is, 

A„ = {neZ: 16n 2 + 16n + 3}, for A ++ _+) and A+ +). Therefore E = E n = 

{neZ: -An 2 - An}. 

The corresponding solutions for A n are 

A+++-) $i.„ = PlnU^lfl, $1,0 = ^~ 1)5 f n = 1, $1,0 = 2 - J, 
A+— + ) $2,„ = Q 2 «/r l $2,0, $2,0 = -T^= fn = 1 ^2,0 = ^=5Tf- 

These two solutions are equivalent to the same solution of the original Schrodingcr 
equation and corresponds to the well known supcrsymmctric quantum mechanics 
approach to this Poschl- Teller potential, [24j[25]. Furthermore, for all A E A„, 

DGal(L A /C(.T)) = G [4] , DGal(2 A /7?) = DGal(i A /A") = e, 

dime £ (H - A) = 2, dim c £{H - A) = dim c £{H - A) = 4. 

Searching Potentials From Parameterized Differential Equations. 

The main object to search new potentials using d z is the family of differential 
equations presented by Darboux in [28] , see section 11.2.21 and equation (|1.25p , 
which can be written in the form 



o2 p , (3z+2)(3z-2) a ft . n(n+8) p _ n 

U z r n T 2 ( 2 + i)( 2 _i) "z-f>i T ( l+ l)(l_,) r n — u > 

«2 p , 3z 2 +2 a p I n(n+2) p _ n 

"I" z(z + l)(2-l) Z ™ + (jH-1)(1-») " _ 

a2p | 6z 2 -3x-4 ^ p | nz(rt+5)+6 g _ n 

"z^" "I" Z(H-1)(Z-1) * ™ Z(^ + 1)(1-Z) " ~~ U ' 

Q2p i 3z-2 a p | nz(n-l)-12 9 _ n 

U z r " "I" z(z + l)(l-z) az ' r ' 1 ^ 4z(z+l)(l-z)- r " ~~ U ' 

«2p i 6z 2 +3z-4 a p i z(n 2 +5)-6 5 _ n 

U z r n "I" 2 ( 2 + i)( 2 _i) c 'z^n "I" 2 ( 2+ l)(l_ 2 )- r " — u ' 

o2p i 3z+2 pt p i nz(n-l)+3 g _ n 

+ 2 ( 2 + l)( 2 _l) C 'z^n + 2 ( 2+1) ( 1 _ z) n i - U, 

u z r n "I- 2 ( z+ i)( z _i) y z r n T ( 2 + i)(i_ 2 ) r n — U, 

^z^n T 2 ( 2 + i)(i_ 2 ) c 'z- r n T ( 2 + i) ( 2 _ 1) — u - 



2 z y + Pd z y + (Q - XR)y = 0, P,Q,Re K. 



(2.13) 
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We recall that some Riemann's differential equations, presented in section [1.1. 41 
corresponds to this kind. 

When we have a differential equation in the form (|2. 13|) , wc reduce it to put 
it in the form of the reduced algebrized Schrodinger equation H$ = A<£>, check- 
ing that Card(A) > 1. Thus, starting with the potential V and arriving to the 
potential V we obtain the Schrodinger equation = A^. This methodology 
(heuristic) is detailed below. 



1. Reduce a differential equation of the form (|2. 13[) and put it in the form 
H$ = A$, checking that Card(A) > 1 and to avoid triviality, a must be a 
non-constant function. 

2. Write W = £ 9* (In a) and obtain V(z) = a(V - d z W - W 2 ). 

3. Solve the differential equation (p x z) 2 = a, write z = z(x), V(x) = V(z(x)). 

To illustrate this method, we present the following examples. 
Bessel Potentials 



• (From Darboux transformations over V = 0) In the differential equation 



we see that A = —n(n + 1) and a = z 2 . Applying the method, we obtain 
V = /i we obtain V(z) = fiz 2 + j and z = z(x) = e ±x . Thus, we have 
obtained the potentials V(x) = V(z(x)) = [ie ±2x + j (compare with 
§6.9]). 

(From Bessel differential equation) The equation 

1 z 2 — n 2 i 
d 2 y + -d z y + —y = 0, n e - + Z, 

z z z 2 

is transformed to the reduced equation 

„ (n 2 4z 2 + l\ 

We can see that A = — n 2 , a = z 2 , obtaining V = — z 2 — |, V = -z a -jZ 2 + j 
and z = z(x) = e ±x . Thus, we have obtained the potential V(x) = V(z(x)) = 
-e ±ix - \e ±2x + \ (compare with [5BJ §6.9]). 
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We remark that the previous examples give us potentials related with the 
Morse potential, due to their solutions are given in term of Bessel functions. 

We can apply this method to equations such as Whittakcr, Hypcrgcometric 
and in particular, differential equations involving orthogonal polynomials (com- 
pare with [H §5]). 
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Final Remark 



The aim of this work is to give, in contemporary terms, a formalization of original 
ideas and intuitions given by G. Darboux, E. Witten and L. E. Gendenshtei'n in 
the context of the Galois theory of linear differential equations. We found the 
following facts. 

• The superpotential is an algebraic solution of the Riccati equation associated 
with a potential, defined over a differential field. 

• Darboux transformation was interpreted as an isogaloisian transformation, 
allowing to obtain isomorphisms between their cigcnrings. 

• We introduced in a general way the Hamiltonian algcbrization method, which 
in particular allow to apply algorithmic tools such as Kovacic's algorithm 
to obtain the solutions, differential Galois groups and Eigcnrings of second 
order linear differential equations. We applied successfully this algcbrization 
procedure to solve problems in Supersymmetric quantum mechanics. 

• We can construct algebraically solvable and non-trivial algebraically quasi- 
solvable potentials in the following ways. 

1. Giving the potential in where for A = Ao the Schrodingcr equation is 
integrable. After we put A ^ Ao checking that the Schrodinger equation 
is integrable for more than one value of the parameter A. 

2. Giving a superpotential to obtain the potential and after we check if 
the Schrodingcr equation is integrable for more than one value of the 
parameter A. 

3. Since parameterized second order linear differential equations applying 
an inverse process in the Hamiltonian algcbrization method. In partic- 
ular, we can use algebraically solvable and algebraically quasi-solvable 
potentials, special functions with parameters (in particular with poly- 
nomial solutions). 

This thesis is a starting point to analyze quantum theories through Ga- 
loisian theories. Therefore open questions and future work arise in a natural way: 
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Final Remark 



supersymmetric quantum mechanics with dimension greater than 2, relationship 
between algebraic and analytic spectrums, etc. 

As a conclusion, as happen in other areas of the field of differential equations, 
in view of the many families of examples studied along this thesis, we can conclude 
that the differential Galois theory is a natural framework in which some aspects 
of supersymmetric quantum mechanics may appear more clearly. 
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